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ABSTRACT 

We adapt the Jain-Taylor (2003) shear-ratio geometric lensing method to measure the 
dark energy equation of state, w = p v /p v and its time derivative from dark matter 
^ ' haloes in cosmologies with arbitrary spatial curvature. The full shear-ratio covariancc 

matrix is calculated for lensed sources, including the intervening large-scale structure 
and photometric redshift errors as additional sources of noise, and a maximum like- 
lihood method for applying the test is presented. Decomposing the lensing matter 
distribution into dark matter haloes we calculate the parameter covariance matrix 
for an arbitrary experiment. Combining with the expected results from the CMB we 
^vq \ design an optimal survey for probing dark energy. This shows that a targeted survey 

■ imaging 60 of the largest clusters in a hemisphere with 5-band optical photometric 
^3 ' redshifts to a median galaxy depth of z m = 0.9 could measure wq = w(z = 0) to a 

marginal 1-cr error of Awq = 0.5. We marginalize over all other parameters includ- 

■ ing w a , where the equation of state is parameterized in terms of scale factor a as 

w(a) = Wo + w a (l — a). For higher accuracy a large-scale photometric redshift survey 

^2 , is required, where the largest gain in signal arises from the numerous « 10 14 M Q haloes 

corresponding to medium-sized galaxy clusters. Combined with the expected Planck 
Surveyor results, such a near-future 5-band survey covering 10,000 square degrees to 
f~| . z m = 0.7 could measure wo to Awq = 0.075 and Aw a = 0.33. A stronger combined 

Qh constraint is put on w measured at the pivot redshift z p = 0.27 of Aw(z p ) = 0.0298. 

We compare and combine the geometric test with the cosmological and dark energy 
^ ■ parameters measured from planned Baryon Acoustic Oscillation (BAO) and super- 

^ ' nova Type la experiments, and find that the geometric test results combine with a 

C^J . significant reduction in errors due to different degeneracies. A combination of geomet- 

ric lensing, CMB and BAO experiments could achieve Awq = 0.047 and Aw a = 0.111 
with a pivot redshift constraint of Aw(z p ) = 0.020 at z p = 0.62. Simple relations are 
presented that show how our lensing results can be scaled to other telescope classes 
',— ] ■ and survey parameters. 
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1 INTRODUCTION 

Over the last decade, gravitational lensing has emerged as 
the simplest and most direct way to probe the distribution 
of matter in the Universe (Bartelmann & Schneider 2001, 
Refregier 2003). More recently it has become apparent that 
it can also be used as a probe of the mysterious, negative- 
pressure "dark energy" component of the Universe which 
gives rise to the observed acceleration of the expansion of 
the Universe (Hu & Tegmark,1999; Huterer, 2002; Jain & 
Taylor, 2003; Hu, 2003; Takada & Jain, 2003; Song & Knox, 
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2004; Ishak, 2005; Ma, Hu & Huterer, 2006; Heavens et al., 
2006) 

The dark energy exerts its influence by its effect on the 
expansion history of the Universe. If the current expansion 
of the Universe is accelerating, the Universe must be older 
than if it was decelerating, since the expansion was slower 
in the past. This changes the distance traveled by a photon, 
r(z), for a given expansion factor of the Universe, a(z) = 
1/(1 + z), as photons have had more time to travel further 
than in the decelerating case. The accelerated expansion will 
also slow the rate of growth of matter perturbations. The 
simplest phenomenological model of the dark energy can be 
constructed by simply parameterizing the equation of state 
of the vacuum, 
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p v = wp v , (1) 

where p v is the dark-energy /vacuum-pressure and p v is its 
energy-density, and w = w(a) may vary with scale factor. 

Gravitational lensing depends upon both the geometry 
of the Universe, via the observer-lens-source distances, and 
on the growth of structure which will lens distant galaxies, 
and so lensing probes both effects. Gravitational lensing is 
an integral effect and so for a given line of sight these effects 
are degenerate with each other and other parameters. In or- 
der to disentangle the effects of the dark energy we require 
redshift information for the source images. It has already 
been shown that such information can be used to reconstruct 
the 3-D distribution of dark matter (Taylor, 2001; Taylor et 
al., 2004). For large-scale imaging surveys, the most practi- 
cal way to get redshifts for each image is from multi-band 
photometric redshift surveys. The COMBO-17 imaging and 
photometric survey (Wolf et al., 2003) has already shown 
the power of combining lensing with photometric redshifts 
(Brown et al., 2003; Taylor et al., 2004; Gray et al., 2004; 
Bacon et al., 2005; Semboloni et al., 2006). 

The parameters of the dark energy can be extracted 
from weak gravitational shear measurements by taking cor- 
relations of galaxy ellipticities at different redshifts (e.g. Ba- 
con et al., 2005; Hu, 2003; Heavens, 2003; Heavens et al., 
2006; Semboloni, 2006), where the expansion history enters 
both the lens geometry and the dark matter evolution rate. 
Jain & Taylor (2003) proposed an alternative approach, tak- 
ing the ratio of the galaxy-shear correlation functions at dif- 
ferent redshifts. In this case the mass of the lens dropped out 
leaving behind a purely geometric quantity useful for mea- 
suring cosmological parameters. This had the advantages of 
allowing the analysis to extend into the nonlinear clustering 
regimes where modeling the nonlinear matter power spec- 
trum can be inaccurate, and where the shear signal will also 
be stronger, i.e. in the vicinity of galaxy clusters. In ad- 
dition, as this relies upon the correlation between galaxies 
and shear many systematic effects will be averaged over, as 
in galaxy-galaxy lensing. Following this a number of papers 
have suggested variations on this theme (Bernstein & Jain, 
2003; Hu & Jain, 2004; Zhang, Hui & Stebbins, 2005). 

Geometric tests of dark energy not only complement 
other methods based on the clustering of matter, but di- 
rectly probe the global evolution of the Universe via the 
redshift-distance relation, r(z). Other methods measure the 
combined effect of the growth rate of perturbations and the 
global geometry. Comparison of the two can be used to test 
the Einstein-Hilbert action, and extensions and modifica- 
tions of General Relativity such as extra dimensions. 

While the main focus of the Jain- Taylor (2003) paper 
was a statistic given by the ratio of galaxy-shear correla- 
tions (or equivalently power spectra), they illustrated their 
method with the analysis of a single cluster. In this paper 
we develop this idea further and focus on applying the geo- 
metric test behind individual galaxy clusters. The main dif- 
ference between this and the original Jain- Taylor approach 
is that we do not need to first generate galaxy-shear cross- 
correlation functions, or cross-power spectra, which require 
large data-sets. Rather the ratios used are just of the shears 
behind a given cluster at fixed redshifts. This allows the test 
to be applied to noisy data, since we do not need to estimate 
correlation functions before applying the ratio test. This is 



similar to the approach of Bernstein & Jain (2004) , who con- 
sidered a "template matching" approach, cross-correlating a 
foreground galaxy template with the background shear pat- 
tern. Our approach is different in that we use the galaxies to 
identify the positions of lensing haloes, and then take shear 
ratios. In doing so we focus on the dark matter haloes gener- 
ating the signal, allowing a halo decomposition of the matter 
distribution, and ask how to maximize the signal. The price 
we pay for this approach is that we become susceptible to a 
sampling variance due to lensing by other large-scale struc- 
ture along the line of sight, which we can beat down using 
multiple lines of sight. In addition we generalize our methods 
to non-flat cosmological models. 

Zhang, Hui & Stebbins (2005) have proposed a different 
geometric method, which allows them to extend the corre- 
lation/power spectrum method to galaxy-galaxy and shear- 
shear correlations as well as galaxy-shear cross-correlations. 
They also point out some inaccuracies with the analysis of 
Jain & Taylor (2003) and Bernstein & Jain (2003), which 
we correct here. 

In the next Section we lay out the basic lensing equa- 
tions we will need. In Section 3 we derive the statistical prop- 
erties of the shear ratios, and write down a likelihood func- 
tion for measuring the dependent cosmological parameters, 
we then estimate the Fisher matrix and parameter covari- 
ance matrix for the dark energy. In Section 4 we outline the 
survey design formalism, using the dark matter halo model 
for the distribution of galaxy clusters and group haloes, out- 
lining a realistic photometric redshift analysis and discuss 
bias and intrinsic ellipticity issues. In Section 5 we discuss 
survey strategies, considering targeted, wide-field and area 
limited designs. Using the parameter covariance matrix and 
a model of photometric redshifts we optimize a weak lensing 
photometric redshift survey for measuring dark energy pa- 
rameters from cluster lensing in Section 6. We forecast the 
expected accuracy of cosmological parameters in Section 7 
and compare and combine with other methods. In Section 8 
we discuss the required control of systematic effects, and we 
present our conclusions in Section 9. We begin by introduc- 
ing the necessary cosmological and weak lensing concepts. 



2 THE DARK ENERGY SHEAR-RATIO 
GEOMETRIC METHOD 

2.1 Background Cosmology 

We start with the metric 

dr 2 = (1 + 2$)dt 2 - a 2 (t)(l - 2$)[dr 2 + St{r)d^ 2 ] (2) 

where r is the invariant proper time, t is the cosmic time, $ 
is the Newtonian gravitational potential, a(t) = (1 + z) _1 is 
the scale factor, r(z) is the comoving distance given by 

and 

{r sin(r/r ) (k = 1) 

r (k = 0) (4) 

r sinh(r/r ) (fc = -1), 

is the angular distance, where ro = 1/ y^hi\Ho is the radius 
of curvature of the Universe and Ho is the current value of 
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the Hubble parameter. The time-variation of the Hubble 
parameter with the cosmic scale factor, H(a), is given by 

1/2 

,(5) 



H(a) 



Q. m a 3 + Q.k cl 2 + Q v e 



-3 J* dlna'[l + w(a')] 



which is a function of the vacuum equation of state, w(a), 
and the present-day density parameters; the matter density, 
Qm, vacuum density, Q v , and the energy-density associated 
with the curvature, Q.k = 1 — f2 m — fi„. A useful expansion of 
the time-dependence of the equation of state in terms of the 
expansion parameter, a, is (Chevallier and Polarski, 2001; 
Linder, 2002) 



w(a) = wo + w a (l — a), 



(6) 



which evolves from w = wo + w a at high redshift to w = wo 
at low redshift, with the transition around z = 1. In this 
case the time-dependence of the Hubble parameter is given 
by 



H\ a ) rn -3 , n -2 , n - 
— — = [\lm,a +\Ikol +il v a 

no 



2.2 Weak Shear 



A galaxy cluster at a redshift zi will induce a shear pat- 
tern on the background galaxies, which can be expressed in 
complex notation as 



7 (0) = 7 i(0)+i72(0), 



(8) 



where 71 and 72 are orthogonal components of the shear field 
at an angle of 6. Around lensing clusters it is convenient to 
use the shear tangential around the cluster centre. This can 
be projected out from the total shear by 



7« = -[71 cos(2<^) + 72 sin(2<p)], 



(9) 



where ip is an azimuthal angle around the centre of the clus- 
ter. 

The amplitude of the induced tangential shear distor- 
tion behind a cluster at redshift z; will grow with redshift 

as 



•Jt(z) = It, 



S k [r{z) - r(zi)] 



zi < z 



(10) 



Zl < Zi < Zj. (11) 



Sk[r(z)] ' 

where 74,00 is the tangential shear induced on a galaxy at 
infinite redshift. If we take the ratio of the shear values at 
two different background redshifts, Zi and Zj, (Jain & Taylor, 
2003) then 

„ _ 7m = S k [r(z j )]S k [r(z i ) - r(zi)] 
lJ Jt,j S k [r(zi)]S k [r{z ] ) - r(zi)]' 

This is the key equation describing the geometric method. 
In the last term the mass and structure of the cluster has 
dropped out. 

In the real Universe galaxy clusters are not isolated, 
and additional large-scale structure along the line of sight 
between the lens and the background source galaxies will 
contribute to the observed shear in both backgrounds. If we 
assume that the large-scale structure is uncorrelated with 
the cluster then this effect will average out over independent 
clusters. Defining 



D, 



7M 
It j 



(12) 
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Figure 1. The response of the shear ratio, R, to each of the four 
cosmological parameters, f2 m , Q v (where = 1 — Q m — Q v ), wo 
and w a , as a function of source redshift. We have rescaled the line 
for w a by a factor l/w a to make it finite, so that dlnR/d w a is 
plotted. Here we have set Z2 = 2z\. The assumed fiducial model 
is H m = 0.27, n = 0.73, w = -1.0 and w a = 0.0. 



as the observed ratio of the tangential shear between two 
redshifts for a given cluster we find on average that 



(Dij) — Rij 



g fc [r(z J )]g fc [r(zQ-r(z i )] 
Sklrizi^SkHzj) - r(zi)] : 



zi < Zi < Zj. (13) 



2.3 Response of Shear Ratios to Cosmological 
Parameters 

The intrinsic sensitivity of the shear ratio, Rij, to a cosmo- 
logical parameter, 9, can be estimated from its logarithmic 
response, 

A9 



ARij _ ( 9 In R^ 



Ru 



(14) 



Figure 1 shows the response of the shear ratio, Rij, to each 
of the cosmological parameters that fix the geometry of the 
Universe for a lens at zi — 0.2, and with backgrounds at 
zi — z and zi = 2z. From this we can see that the response 
of the shear ratio to cosmological parameters is weak, scaling 
roughly as 



D 1 i-O.CWoO.OloO'OtW -0.0()liu a 

R oc \w \ il v S2 m e 



(15) 



for sources at z\ = 1 and Z2 = 2. This weak dependence calls 
for high accuracy in the shear measurements. We discuss the 
control of systematics in Section 7. 

The similarity of the responses of the shear ratio to 
different cosmological parameters in Figure 1 also indicates 
their strong degeneracies. We can expect that wo will be 
correlated with Q v and Q m , whilst the Q v -Q m and wo- 
w a combinations will be anti-correlated with each other. 
The similarity of the responses of Wo and fi m suggest 
these parameters will be highly degenerate, while the dif- 
ferences between wo and il v at low redshift suggest these 
should be less correlated. The response of R to wq peaks at 
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Figure 2. The tangential shear profile as a function of wq for 
a lens at z = 0.2 normalized relative to 7(2 = 00), showing the 
effect of any shape changes. The lines are, from lowest to highest 
are for wq = —1.5, —1.0, —0.5. 
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Figure 3. The tangential shear profile as a function of w a for 
a lens at z = 0.2 normalized relative to 7(2 = 00), showing the 
effect of any shape changes. The lines are from lowest to highest 
are for w a = -0.5,0.0,0.5. 



around z — 0.8, when the dark energy begins to dominate 
the energy-density of the Universe. Interestingly, although 
weak, the geometric shear ratio method is most sensitive to 
Wo. As we parameterize the dark energy equation of state as 
w(z) = wo +w a [z/(l + z)] a change in wo affects the ampli- 
tude of w(z) at all redshifts and hence affects the shape of 
the tangential shear as a function of redshift at all redshifts. 
At low redshift w a only changes the slope of w(z) and its 
amplitude at higher redshift where the effect of dark energy 
is less significant. 

Since we take shear ratios, we are only sensitive to 
changes in the shear-redshift relation. Figure 2 and 3 show 
the shear as a function of z normalized to unity at z = 00 



using equation (10) for different wo and w a - This shows the 
effect that changes in wo and w a have on the shape of the 
tangential shear. Varying both wo and w a by 0.5, we see 
that wo has a much larger effect than w a on the shape of 
the shear as a function of redshift. 



3 STATISTICAL PROPERTIES 

In this Section we present a maximum likelihood approach 
to measuring the geometry of the Universe from shear ratios 
around individual clusters and galaxy groups. In this analy- 
sis we shall consider shot-noise, from galaxy discreteness and 
intrinsic galaxy ellipticities, the effect of lensing by large- 
scale structure between the lens and the background source 
galaxies which will act as an addition source of correlated 
and uncorrelated noise, and photometric redshift errors. 

3.1 Likelihood Analysis 

Compressing the notation for a pair of background galaxies 
as n = and v = (m,n) we can write the covariance 

matrix for shear ratios as 

= (ARvARp). (16) 

The log-likelihood function for the four cosmological param- 
eters estimated from a single cluster is then 

— 21nL c (flu,f2 m , wo, w a \D) = 

Y J i.R^-D li )\C^]-\R v -D u ). (17) 

Here we have further assumed that the scatter between R v 
and D v , due to shot-noise, photometric redshift errors and 
cosmic shear from large-scale structure, is Gaussian dis- 
tributed. We do not need to assume the lensing signal from 
the clusters itself is Gaussian. 

If our survey contains multiple independent clusters, the 
total log-likelihood is just the sum of the log-likelihoods for 
the individual clusters; 

lnL T OT(fi 1 >,fim,wo,Wa) = In L c (Q v , Q. m , w , w a \ D) , (18) 

c=l 

where N c \ is the number of independent clusters in the sur- 
vey. 

3.2 The Covariance of R 

The covariance matrix of shear ratios, C^f 1 , is given by 

(ARij AR mn } — (ARijAR mn } Bn + (ARij AR mn }i ssy (19) 

which can be decomposed into a shot-noise term due to the 
intrinsic dispersion in galaxy ellipticities, and a term due to 
cosmic shear induced by the intervening large-scale structure 
between the lens and the two background sources. 

To avoid double counting and taking ratios of the same 
redshift bins, the indices in equation (19) are restricted to 
i < j and m < n. Finally, there is a remaining degeneracy 
between the shear ratios, since the ratios between any three 
galaxies at redshift Zi, Zj and z^ obey the relation 

/>'., //.,•.///.... (20) 
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This reduces the total number of permutations of usable 
pairs of galaxies to (N g — 1), where N g is the total number 
of galaxies. In practice we will bin data, in which case this 
also applies to bins. 



3.2.1 Shot-Noise Covariance 

The first term in equation (19) is due to shot-noise, arising 
from the discrete nature of galaxies and the intrinsic disper- 
sion in galaxy ellipticities; 

(ARij sn 



7:' ' V 



K 



+ ( ^ 

1 Jj 



iOjn Ojm J t 



where 5^ is the Kronecker delta-function, and 



(21) 



(22) 



is the fractional variance in the tangential shear due to 
the intrinsic dispersion in background galaxy ellipticity per 
mode, <r e , and 7, = y(zi) is the expected tangential shear 
signal from the cluster for a background galaxy at redshift 
Zi. Here we shall use <r e = 0.3 per mode. 

There is a subtlety in determining the distribution of 
the fractional variance for the ratio of two ellipticity mea- 
surements. If we assume that the observed ellipticities have 
zero mean, and that the distribution of intrinsic galaxy el- 
lipticities is Gaussian, the resulting distribution of the ratio 
of ellipticities has a Cauchy/Lorentzian distribution, and so 
an infinite variance. Around a lensing cluster, the mean el- 
lipticity of the background galaxies is non-zero, and if we 
assume that the mean signal is always greater than the dis- 
persion in the mean due to intrinsic galaxy ellipticities the 
variance is finite and we can assume Gaussian errors. This is 
certainly the case when we average the shear both tangen- 
tially around a cluster and in redshift bins. Hence, instead 
of working with individual galaxies we shall consider binned 
galaxies, where the fractional variance in the shear is 



7» 



(23) 



per redshift bin, where Ni is the number of galaxies in the 
i th redshift bin. From hereon, the indices i, j will refer to 
bin number, rather than individual galaxies. 



3.2.2 Photometric Redshift Errors 

In current and future weak lensing surveys photometric red- 
shifts will also be available as an estimate of galaxy dis- 
tances (see, e.g., the COMBO-17 photometric redshift sur- 
veys, Wolf et al, 2001; CFHTLS, Semboloni et al., 2006). 
Here we characterize the effect of photometric redshift un- 
certainty on shear ratios. 

The effect of errors on the photometric estimates of 
galaxy redshifts is to dilute the shear signal in each red- 
shift bin by randomly moving galaxies in and out of any 
particular bin. If we assume that the distribution of redshift 
errors is a Gaussian with width o z (z g ) which depends on the 
true redshift of the galaxy, z g , and has a bias in the mean 
of the distribution Zbi as , then 



p(z\z g ,o z ) = 



-(z-z g + z bias ) /2<T z (z g ) 



(24) 



We shall take Zbias = for all experiments, but its effect on 
the marginal error of w(z) will be discussed in Section 4.3, 
where we a lso discuss the e ffect of a change in the variance 
a z {z) — ► y/ (t'1(z) + Aerf (z). We discuss the specific form for 
<J z {z) for photometric redshift surveys in Section 4.2. 
The expected shear in a redshift bin is given by 



(7t,i) = It,, 



dzn(z) Sk[r i Z ]„ f J^ Zl)] iM*-*|g.(st)],(25) 



Sk[r(z)] 



where (e.g., Ma et al., 2005) 



Paz[z\(t z 



erf 
erf 



Z + Zbias + Az/2 
Z + Zbias — Az/2 

Vlo-z 



(26) 



is the part of the redshift error distribution which lies in 
a redshift bin of width Az centred on z, and erf(a;) is the 
error function. The estimated shear is weighted by the num- 
ber of galaxies scattered from one redshift to another, given 
by the galaxy redshift distribution, n(z). Equation (25) is 
normalized so that 



dzn(z)P&z[zi — z\o z (z) 



for each redshift slice at Zi 



3.2.3 Cosmic Shear Covariance 



(27) 



The second term in equation (19), due to the cosmic tan- 
gential shear induced by large-scale structure between the 
lens and the source planes, is given by 



(ARij A7? mn )i ss _ Ci t 



im 



Rij Rmn 



+ 



s~t~it "ft f^"ft "it 



Jt,ijt,m 7*,j7t,n 7t,i7t,n 7t,i7t,m 



+ 



rj~tt "ft 

2,i,min(j,n) S K 



Qltlt 

„ K 2,max(i,m) ,j <-K . , 
"im i — °jn\ Z °) 



with the same restriction on indices as for the shot-noise 
term. Ci and Ci are defined in equations (29), (30) and 
(31). The first four terms in equation (28) are due to the 
correlated distortions induced on both background galaxy 
images by matter lying in front of the nearest source plane. 
The last two terms arise from matter lying between the back- 
ground source planes and can be regarded as an extra noise 
term on the ellipticities of the furthest background source 
galaxies. The covariance of the induced tangential shear for 
background galaxies at redshifts Zi and Zj due to large-scale 
structure between the observer and the background source 
galaxies and averaged over an aperture of radius 6 is 

2[1-J O {£0)] ,h(£8) \ 2 

P6 2 £6 J ' 1 ' 

where a = (1, 2). This is derived in Appendix A. Here J n is 
the n th order Bessel function, and the angular shear-shear 
power spectrum for the two source galaxies is 



C. 



«(«)= / -r c ir 



and 



r i <r j 



dr P s [£/S k (r),r]W[r, n]W[r, rj ],(30) 
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drP s [£/S k (r),r]W 2 [r,i 



(31) 



for sources at redshifts Zi and Zj, and r% = r(zi). We have 
used a nonlinear matter power spectrum, Ps(k, r), with a 
ACDM model with concordance parameter values, Q m = 
0.27, f2„ = 0.73, ft = 0.71, using the functional form of 
Eisenstein & Hu (1999) for the linear power spectrum. The 
linear power spectrum is mapped to the nonlinear regime 
using the fitting functions of Smith et al. (2003). The lensing 
weighting function in equations (30) and (31) is given by 



W{r,n] 



S k (r 



Sk(n)a(r) 



(32) 



In the case of binned data with photometric redshift errors, 
this becomes 

/>oo 

W(r,n)= dzn{z)P Az [ Zl - z\a z (zi)]W(r,r(z)). (33) 

J z, 

These integrals are normalized as in equation (27). 



3.3 Parameter Covariances 

This likelihood analysis can be used to extract cosmological 
parameter error estimations. The parameter covariance ma- 
trix can be calculated from the inverse of the Fisher matrix. 



where 

d 2 InL 



86r 



(34) 



(35) 



is the Fisher matrix and 9 = (Q v ,Q m ,wo,w a ) is a vector 
containing our cosmological parameters (see Tegmark, Tay- 
lor & Heavens, 1997, for an introductory review). 

For a Gaussian Likelihood function with parameters in 
the mean, such as equation (17) the Fisher matrix reduces 
to (e.g Tegmark, Taylor & Heavens, 1997) 

F ^ = \ J2 (OtRn[C™]~ 1 d i R* + dMC™]- 1 ^), (36) 



where di denotes differentiation in parameter space, and 
the summation in n and v denotes summing over all non- 
degenerate source configurations (see Section 3.2). The 
marginalized error on the parameters is given by 



<A0?> r 



= \F- 



while the conditional error is 
<A0 2 ) cond = 1 1 Fa < [F-%. 



(37) 



(38) 



Throughout we shall quote marginalized errors. Results on 
parameter accuracies are presented in Section 7. 



4 SURVEY DESIGN FORMALISM 

To understand the contribution to the geometric test signal, 
we use a halo decomposition of the matter density distri- 
bution (Peacock & Smith, 2000; Smith et al., 2003; Seljak, 
2000). The full signal then comes from integrating over all 
halo masses, lens redshifts and background sources, but with 
the halo decomposition we can extract information about 



which halo mass range contributes most to the signal. This 
will help to determine optimal survey strategies. 

In this Section we shall also discuss a more detailed 
model for photometric redshift errors, based on studies of 
photometric redshift accuracies from the COMBO-17 sur- 
vey (Wolf et al., 2003), and the limits of ground-based mea- 
surements of galaxy ellipticities. These elements are then 
factored into the optimization of a weak lensing survey in 
Sections 5 and 6. 



4.1 Halo Decomposition of the Matter Density 
Field 

So far we have only considered the shear signal from a single 
cluster. In practice we would want to sum over many galaxy 
clusters in a weak lensing survey. In this case we need a 
model for the abundance of clusters as a function of mass 
and redshift, M{M, z). To apply this we must first find the 
relation between mass and shear. For simplicity we shall use 
the singular isothermal sphere model. 

The mean shear signal inside a circular aperture of an- 
gular radius 6 for a singular isothermal sphere, and a source 
with virial mass M at infinity is 



%, ao {<0,M) = 



AM) 



where 



,(M) 



- ( M \ 



(39) 



(40) 



is the Einstein radius for a source at infinity. In the last 
expression we have made use of the constant virial velocity 
of the singular isothermal sphere, 



3GM 
2r v ' 



(41) 



where the virial mass, M, is the mass enclosed by the virial 
radius, r v ; 



M = -?-rlp m 5 v . 



(42) 



Here p m is the mean mass-density of the Universe, S v = 
340 (Eke et al. 1996) is the virial overdensity for a ACDM 
Universe, and 



M 



(43) 



is a characteristic mass. Equation (42) defines the virial ra- 
dius, r v , in terms of the virial mass, M, and is given by 



r v = 0.293 



/ M 



\ 10 13 A/ G 



1/3 



(O m /i 2 ) 1/3 Mpc. 



(44) 



Substituting this into the expression for the velocity dis- 
persion, g v , in equation (41) and then into the expression 
for 0oo (equation 40) we find the shear signal scales as 
7t oc A/ 2//3 . As more massive clusters are larger their sur- 
face mass-density, and hence mean shear, scales more slowly 
than in proportion to the mass, as would be expected for 
fixed sized haloes. 

The shot-noise term in the shear covariance matrix for 
sources at infinity is 
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A 7 2 

7t,oc 



(45) 



where n z (> z{) is the surface density of galaxies that lie at 
a redshift greater than the lens at zi. The angular radius, 
6, drops out of this expression since the signal-to-noise ratio 
for the mean shear of an isothermal sphere is a constant for 
a uniformly distribution of background sources. Hence the 
Fisher matrix will scale as Fij oc 6%^ oc 

The Fisher matrix for the halo geometric test, integrat- 
ing over lensing cluster mass, M, and lens redshift, zi, is 



p. 



J dzi 



dU 



I M \ 4/3 



fi J (Mo,«,)(46) 



where have factored out the mass-dependency of the Fisher 
matrix and set Q x = 1 in Fij (Mo, Zi) for a single halo. In this 
expression M(M,zi) is the number density of clusters per 
[ft. _1 Mpc] 3 with mass M at redshift zi. The lower mass cut- 
off in the integral over mass, M_ (zi), is set by the condition 
that a cluster shear must be measurable with a signal-to- 
noise of 



7 

which sets 
M-(zi) =M 



-i 3/2 



a e fi 



\/n z (> zi)n 



(47) 



(48) 



We shall assume fi = 1 from hereon. Note that although we 
have a low signal-to-noise threshold for measuring the shear 
signal from a given halo, we assume that the detection of a 
halo is based on the detection of galaxies in the halo, and 
therefore has a high signal-to-noise. 

The halo number density is a function of mass, M, and 
redshift, z, given by 



J\f(M, 



M 



f(M,z) 



(49) 



where fraction of matter, f(M,z), in haloes of mass M at 
redshift z, can be written in the universal Sheth-Tormen 
form (Sheth & Tormen, 1999), 

-0.3 



1/2 -v/2\/2 



vf(y) = B (l + v/V2) (f/V2) 
where B is a constant of normalization so that 



dvf{v) = l, 
to 

and 

<5 2 

v ■ 



a 2 {M,z) 



(50) 



(51) 



(52) 



The form of equation (50) finds justification from the el- 
lipsoidal collapse model of haloes (Sheth et al, 2001). The 
collapse threshold for linear matter overdensities, 5 = <5p/p, 
is 8 C = 1.686. The variance of overdensities in spheres of 
radius R(M) is 

j3 , 



a(M,z) 



where 



d A k 
(2tt)3 



Pun(k,z)j^[kR(M)], 



Pnu{k,z) = D 2 {z)Pnn{k,z = G) 



(53) 



(54) 




Redshift 

Figure 4. The cumulative number count of dark matter haloes 
per square degree, J\f(> M, < z), as a function of redshift, for the 
mass range M = 1O 12 M to M = 10 15 M Q . 



is the linear matter power spectrum. This can be split into 
a linear growth factor, 

(l + z')dz' 



D(z) = H(z) 



H(z') 



(55) 



(see Heath, 1977; Carroll, Press & Turner, 1992; Linder, 
2003) where H(z) is given by equation (5), and the present- 
day linear matter power spectrum, Pn n (fc, 2 = 0). This equa- 
tion is only valid for wo — { — 1, —1/3, 0}. We use it here 
as ACDM is our fiducial cosmology. The spherical Bessel 
function, jo{x) = sin(a;)/a; is the transform of a sphere. The 
radius of the sphere, R(M), is the linearized radius of a clus- 
ter halo of mass M, which will collapse down to a nonlinear 
virial radius, r v , and is given by 



R(M) = [1 + z)Sl /s r v (M). 



(56) 



We plot the cumulative number count of dark matter 
haloes per square degree as a function of redshift, for a range 
of halo masses in Figure 4. Typically, for 1O 12 M0 haloes we 
expect 10 3 haloes per square degree, while for 1O 15 M0 haloes 
we expect 10 -2 haloes per square degree. 

In Figure 5 we show the expected cumulative number 
count of dark matter haloes for a range of median redshifts, 
Zm, 7V(> M,z m ) per square degree. The dotted lines rep- 
resent various upper redshift limits with no signal-to-noise 
limit on cluster detection for z m = 0.7. The solid lines are 
for a detection threshold of clusters with signal-to-noise of 
unity for various median redshifts, given by equation (48). 
In assuming a SIS model for the lensing clusters the aver- 
age shear around a cluster may be systematically underes- 
timated. A more reliable model is the Navarro-Frenk- White 
(NFW) profile, although the density profile form would yield 
a more complex relation for the shot noise term than the SIS 
profile. There is, however, an approximate scaling relation 
which relates 7s/s to ^nfw outlined in Wright & Brainerd 
(2000) which, since we take the average tangential shear in 
an apeture, should be adequate. Adopting the techniques 
outlined in Wright & Brainerd (2000), the concentration 
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Figure 5. The cumulative number count of dark matter haloes for 
a range of median redshift distributions, jV(> M, z m ), per square 
degree. The solid lines assume a maximum redshift in the halo 
population of z ma x = 1-5, while the upper dotted line assumes 
£max = 2.5, and lower dotted line z max = 0.5. The cut-offs in 
halo numbers for the z max = 1.5 (solid) lines are for different 
median rcdshifts with a shear signal-to- noise limit \x > 1. 




log 1Q (Mass/M ) 

Figure 6. The ratio of mean shears for an SIS and NFW haloes 
of varying mass for haloes at a redshift of z c = 0.1. The solid 
line is for a ACDM fiducial cosmology, the dashed line is for a 
SUGRA fiducial model and the dot-dashed for a Phantom model, 
see Section 7.5 for details. 



parameter depends on the mass, redshift and fiducial cos- 
mology. We use the concentration parameter from Dolag et 
al. (2004), our scaling from jnfw to jsis depends on the 
dark energy fiducial model, mass and redshift of the cluster. 
We will use this scaling to correct the shear signal expected 
from the halo model. Note that this only affects the noise 
properties since the shear ratio only depends on the redshift- 
distance relation. Figure 6 shows how the scaling depends 
on mass and the fiducial dark energy models (discussed in 
Section 7.5) for clusters at a redshift of z c — 0.1. 



If the haloes are assumed to be randomly distributed 
over the sky, and we take their physical size to be the virial 
radius, the effect of overlapping haloes projected onto the 
sky is negligible. For instance, a M — 10 15 Mq halo has 
a virial radius of r = 0.75Mpc and a number density of 
n ~ 10 -2 per square degree, while a M = 1O 13 M0 halo 
has a virial radius of r = 0.15Mpc and a number density 
of 7i ~ 10 2 per square degree. At z — 0.2, the physical 
distances 0.75Mpc and 0.15Mpc subtend 0.12 degrees and 
0.025 degrees respectively. Hence we shall assume that halo 
overlaps are not important. 

4.2 Photometric Redshift Uncertainty 

In Section 3.2.2 we introduced the effects of including photo- 
metric redshifts on the lensing measurements. Here we detail 
our estimate of the photometric redshift errors. 

The uncertainty on the photometric redshift error on 
an individual galaxy with redshift z and magnitude R for a 
multi-band survey is well fitted by (Wolf et al. 2004); 



(z, R) = A(l + z) [l + lO^*-** 5 ] 



1/2 



(57) 



where A = 0.035, B = 0.8 and R* = 23.0 for galaxies in a 
5-band survey, and A = 0.007, B = 0.8 and R, = 21.6 in a 
17-band, COMBO-17-type survey. This shows that the red- 
shift errors are well constrained at bright magnitudes but 
poorly constrained at faint magnitudes. The first parame- 
ter, A, characterizes the best performance achievable in the 
bright domain, where photon noise is irrelevant and spectral 
resolution limits the redshift estimate. The second parame- 
ter, B, describes how a decrease in photon signal propagates 
into the redshift signal. This should be 0.8 if we consider all 
galaxies, but can be made smaller by filtering out galaxies 
with outlying redshift errors. The final parameter, R„, deter- 
mines the magnitude where we see a sharp rise in the redshift 
error function when we change from the spectral-resolution 
limited regime at bright magnitudes into the regime where 
photon noise drives the redshift noise by a factor of ~ 2.5 per 
magnitude under the assumption of a locally linear transfor- 
mation from colour-space into redshift-space. 

The average redshift error in a bin at redshift z is given 
by averaging over all observable galaxies below a limiting 
absolute magnitude in that bin, 



<J z (z) = 



dMt>(M)a z (z,M) 



I. 



dM *(M) 



(58) 



Here \t(Af) is a sum of Schechter functions <3? rG d and "I'biue 
(see Wolf et al, 2003 for details of the COMBO-17 lumi- 
nosity functions) for a red and blue sample of galaxies. The 
luminosity functions are defined for a colour, c, as 

$ c (M)dM = 0.41nl0^X c Q (M)e" x(M) dM, (59) 

where 

X(M) = io-°- 4 < M -"c*), (60) 
and 

0rcd(^) = (2.0- z) x 10 _3 [/i _1 Mpc] -3 , (61) 

^biue(^) = 3.0 x lO-^/^Mpc]- 3 , (62) 
valid for z < 2, are the characteristic space-densities of 
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Redshift 

Figure 7. Variation of cr z (z) with redshift for a 5-band (upper 
solid line) and a 17-band (lower solid line) photometric redshift 
survey, averaging over galaxy luminosities, for a survey with me- 
dian redshift z m = 0.7 (solid line). Galaxy properties are from 
COMBO-17 and described in the text. Also shown is a standard 5- 
band photometric redshift model with u z (z) = 0.05(1+2) (dashed 
line). 

galaxies. The slope of the luminosity functions are 



«rcd = —0.5, (63) 

abiuo = -1.3, (64) 

and 

M* ed (z) = -20.18 - 1.04z, (65) 

Mu ue (z) = -20.09 - 1.28*, (66) 



are the characteristic magnitudes of reds and blue galaxies 
in the COMBO-17 survey. Mn m (z m ) is the limiting apparent 
magnitude of survey with median redshift z m given by (see 
Brown et al. 2003, and equation 79 in Section 4.5) 

M Vim = 20.8 + z m /0.23 (67) 

for an optical survey, which we then transform to the abso- 
lute limiting magnitude; 

Mum = Mi im - 51og 10 {(1 + z)S k [r{z)]} + K[z). (68) 
The K-correction, K(z), is; 

K(z)=2.5(y-l)]og 10 (l+z), (69) 

where v is the spectral slope of galaxies. We take this to be 
v = +1, making the K-correction zero. Figure 7 shows the 
increase in mean photometric redshift uncertainty for a 5- 
band and a 17-band survey with median redshift z m = 0.7, 
based on the galaxy luminosity functions. As the magnitude 
of a galaxy depends on its redshift, the scaling of the pho- 
tometric redshift noise is more complicated than the simple 
(1 + z) scaling commonly used. Brodwin et al. (2003) find 
<J z (z) — 0.05(1 + 2:) for a 5-band survey, which we plot as the 
dashed line in Figure 7. Our estimate of the redshift error 
for a 5-band survey predicts a higher error for z > 0.7, and a 
lower error for z < 0.7. We have extrapolated these formulae 
to z — 1.5 though this extrapolation may be optimistic as 




i_, 1 , , , , 1 , , , , 1 , , , 1 

0.5 1 1.5 

Redshift 

Figure 8. The effect of photometric redshift errors on the tangen- 
tial shear behind a lensing cluster of mass 10 15 Mq at a redshift of 
0.2, assuming 5 bands. The dashed line is the true shear response, 
while the solid line is the shear with photometric redshift errors, 
using equation (25). 

photometric redshift estimates can increases dramatically at 
z m 1 if IR data is not available. 

For an intermediate 9-band survey we linearly interpo- 
late between the 5-band and 17-band lines, assuming that at 
each redshift the relationship between bands is linear. Over 
all redshifts we find there is no simple linear scaling relation 
with (l + z). However we find approximate fitting formula 
for a 5-band survey, 

ctz(z) w 0.063(0.64 + 2), (70) 
and for a 17-band survey, 

a z (z) w 0.041(0.37 + 2). (71) 

Figure 8 shows the effect of a 5-band photometric redshift 
error, given by equation (25), for a photometric galaxy sur- 
vey parameterized the same as the COMBO-17 survey, with 
median redshift z m — 0.7 and limiting magnitude R — 24, 
on the measured tangential shear distribution behind a 
M = 10 15 M© halo at 2 = 0.2. The main effect is a sup- 
pression of the shear signal at low redshift, where the shear 
is rapidly changing. This is due to the scattering of unlensed 
galaxies in front of the lensing halo into bins just behind the 
halo. 

The photometric redshift error fit from COMBO-17, 
given by equation (57), is per galaxy. In practice the pho- 
tometric redshifts produced by any multi-band analysis will 
also provide an individual redshift error for every galaxy 
which will also depend on redshift and magnitude. In this 
current analysis photometric redshifts are averaged over all 
galaxy types and magnitudes. In practice one would like to 
weight the data optimally to minimize the effect of both 
shear and photometric redshift errors. Given the redshift 
dependence of the shear signal behind a lens, it is likely 
that both errors in the shear signal and photometric red- 
shift errors degrade the measurement of parameter, while at 
redshifts far from the lens, shear errors will dominate. This 
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implies that there is an optimal weighting scheme which is 
a function of galaxy redshift and magnitude for weak shear 
analysis using photometric redshifts. We shall explore this 
elsewhere. 



4.3 Bias in the Photometric Redshifts 

In addition to the uncertainty on photometric redshifts, we 
would also like to know the effect of a bias in the photometric 
redshifts, leading to an off-set in their calibration. We can 
model the effects of this by considering the first-order effect 
of such a bias on the measurable parameters. In Appendix B 
we show that for a Gaussian distributed likelihood function, 
the linear bias in a parameter, which we shall call S9i, due 
to a bias in a fixed model parameter (i.e., one whose value 
we have assumed and is not being measured), which we shall 
call Sipj, is given by (see also Kim et al., 2004) 



(72) 



where F ee is the parameter Fisher matrix and F 6 ^ is a 
pseudo-Fisher matrix of derivatives with respect to parame- 
ters which are assumed fixed (ip) and those to be determined 
(0). 

Assuming there is a possible bias in the mean of the pho- 
tometric redshifts of the survey, Zbias, (see Section 3.2.2) due 
to poor calibration of the photometric redshifts with spec- 
troscopic redshifts, and marginalizing over all other cosmo- 
logical and dark energy parameters, we find that the induced 
bias in wo due to the bias in galaxy redshifts is 



Sw = — Cbias Sz hi , 



(73) 



where Cbias is a constant. If the bias in the mean of the 
photometric redshifts arises from an overall bias in the pho- 
tometric redshift calibration, the calibration error will be 



c( 2 bia 



*(*) 



(74) 



where N spcc is the number of galaxies with a spectroscopic 
redshift. If we set febias = cr(2bias) and a requirement that 
the bias in wq is half of the error, 8wo = 0.5Aiuo, then the 
number of galaxies with spectroscopic redshifts we require 



N spec — 



CbiasCr(«) 



Swo 



(75) 



We have found that Cbias ~ 9.0 for the geometric test. If 
we assume o(z) ~ 0.1 and Awo ~ 0.01 then we require 
TVspoc « 3 x 10 4 . The size of the required spectroscopic red- 
shifts required to calibrate the geometric test suggests that 
a large spectroscopic survey, such as that proposed for the 
Wide-Field Multi-Object Spectrometer (WFMOS; Bassett 
et al., 2005), would be required and combined with a large- 
scale weak lensing survey. 

We have also investigated the effect of an offset in 
the variance of the photometric redshift errors ct z (z) — > 
c| (jz) + Act? (z). We find that this effect is negligible for 
the geometric test, so that the total bias due to photometric 
redshift errors is only dependent on the bias in the offset of 
the mean. However in the pseudo-Fisher analysis the vari- 
ation about the mean of Aa z (z) is ±0.05, and we would 
expect there to be an effect at some level if the variation 



was larger. We explore fully marginalizing over nuisance pa- 
rameters in a full Fisher analysis elsewhere. 

4.4 Limits on the Measurement of Galaxy 
Ellipticity 

4-4-1 Ground-based Ellipticity Measurements 

For ground-based weak lensing observations estimates of 
galaxy ellipticities are limited by atmospheric seeing. The 
angular sizes of typical galaxies in the GOODS fields scale 
with redshift by (Ferguson et al. 2002) 

9 g — 0.8Z 1 arcseconds. (76) 

If S is the typical seeing during weak lensing observations, 
the post-seeing galaxy image will be 



VffJ+el 



(77) 



This will tend to decrease the ellipticity of galaxy images. 
Much effort is put into weak lensing to correct this effect. 
However, once the seeing disc exceeds the galaxy image and 
g <C 6 3 , this correction fails. Typically, galaxy sizes are 
about 6 g — 0.8 arcseconds at redshift z = 1. If the ground- 
based seeing for weak lensing is typically 8 S m 0.7 arcsecs, 
then by a redshift of z = 1.5, the galaxy sizes have dropped 
to 6 — 0.5 arcseconds and galaxy ellipticities cannot be re- 
covered without the use of adaptive optics. 

Another limitation which could potentially come into 
play is when the galaxy image is too faint to properly mea- 
sure the galaxy shape against the sky background. However, 
Bacon et al. (2001) find that the dispersion on the measured 
galaxy ellipticities is very insensitive to the galaxy magni- 
tude, and seems only limited by the detection threshold for 
galaxy detection. For 5-ct detected galaxies, ellipticities can 
be measured down to the limiting magnitude of the survey, 
with ct e = 0.3. 

Given these two results we shall assume that we cannot 
measure redshifts beyond z = 1.5 from the ground due to 
being unable to recover the pre-seeing ellipticity. 

4-4-2 Space-based Ellipticity Measurements 

Rhodes et al. (2003) find no dependence on ellipticity disper- 
sion as a function of redshift for space-based data. Refregier 
et al. (2003) and Massey et al. (2004) find that a e = 0.2 
is a reasonable measure for the ellipticity dispersion for a 
space-based weak lensing survey. They also find a maximum 
redshift bound for space-based surveys can be set at z — 2.0 
corresponding to a deep magnitude cut of R = 29.1. 



4.5 Optical Surveys 

In this Section we outline how to parameterize a weak lens- 
ing and photometric redshift survey, and how these will scale 
for different telescopes. A reasonable way to compare be- 
tween potential survey designs is to consider equal-time ob- 
servations. Hence one can compare dark energy results both 
for a single telescope class, and across telescope classes. The 
time taken for an imaging survey on a given telescope scales 
as (cf equation 88) 



T oc zf n f sky D 2 (fov/l c 



(78) 
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where D is the diameter of the primary mirror of the tele- 
scope and fov is its field of view. We normalize the timescale 
of a survey to the 5-band (g', u, r', i', z') CFHT survey, where 
T = 162 nights for z m = 1.17 (r = 25.9), / sky = 4.25 x 1(T 3 , 
D — 3.6m and fov = 1 square degree. The median redshift 
for an R-band survey is (Brown et al., 2003) 

z m = 0.23{R- 20.6), (79) 

while we find that the projected surface number count den- 
sity of galaxies in the COMBO-17 survey scales with the 
median redshift as 

Ti2(z m ) — 30 z^ 4 galaxies per square arcmin. (80) 

We also need to assume a functional form for the galaxy 
redshift distribution which we take to be 

p(z\z m ) oc z 2 cxp [- (z/z*) 1 ' 5 ] , (81) 

where z* = z m /1.412 and 

/ dzp(z\z m ) = 1. (82) 
Jo 

The space density of galaxies as a function of galaxy redshift, 
z, and survey median redshift, z m , is then 

n,3(z\z m ) = n2{z m )p(z\z m )- (83) 

The 3-D galaxy redshift distribution, n(z) = ns,{z\z m ), is 
used in equations (25) and (33) when calculating the effects 
of photometric redshift errors, for calculating the number 
of galaxies in redshift bin, Ni, for the shot- noise, and for 
finding the cumulative surface density of galaxies above a 
halo redshift, ri2(> z), in equation (45). 

The number of redshift bins used in the background to 
the lenses, Nb, is determined by the photometric redshift 
uncertainty (Section 4.2) by assigning a bin width at par- 
ticular redshift to be the average photometric uncertainty, 
ct z (z), at that redshift. The bins exhaustively fill the avail- 
able redshift range. 



5 SURVEY DESIGN STRATEGY 

Having formulated the basic method for estimating dark en- 
ergy parameters from shear ratios, we now consider the prob- 
lem of what type of survey would be optimal for measuring 
the properties of the dark energy from the shear ratio geo- 
metric test. For instance, should one construct a wide area, 
but shallow, multi-band survey, or a narrow and deep multi- 
band survey with a survey-class telescope, such as the VST 
(Belfiore et al, 2005), the Dark Energy Survey on the CTIO 
(Wester, 2005), darkCAM (Taylor, 2005) or Pan-STARRS 
(Kaiser, 2005)? Or should one instead take snap-shots of 
galaxy clusters with a large but small field-of-view tele- 
scope such as SUBARU, the VLT or the Keck Telescope? 
We shall compare these different strategies by minimizing 
the marginalized uncertainty on wo for fixed-time observa- 
tions, assuming a prior from the expected 14-month Planck 
Surveyor experiment (Lamarre et al., 2003) to lift the main 
degeneracies. 

Broadly we have two observing strategies available to 
us: targeted observations at individual clusters, and a gen- 
eral wide- field survey. In the former, one would use a large 
telescope with small field-of-view to take rapid observations 



of each cluster, while in the latter a large telescope with 
a wide field-of-view would make a general wide-field survey 
from which one would extract haloes. With a halo decompo- 
sition analysis of the matter distribution we can ask where 
most of the signal will come from for a dark energy analysis 
with weak lensing, and so see which strategy would be most 
effective in terms of telescope time. We begin with targeted 
observations of clusters. 



5.1 Targeted Observation Mode 

We shall assume we have a large telescope with a small field- 
of-view which can target pre-selected galaxy clusters from 
a pre-existing galaxy cluster catalogue. The survey would 
start by imaging the largest clusters on the sky, and then 
move on to subsequently smaller haloes. We shall assume 
that the telescope has some fraction of the sky available to 
it, and we shall ignore scheduling issues. 

Figure 9 shows the accuracy on wo, marginalized over 
Q m , fit; and w a , which can be achieved by a targeted sur- 
vey as a function of the number of clusters in decreasing 
mass. We have assumed half the sky (20, 000 square degrees) 
is available, and combined the lensing result with a 4-year 
WMAP prior (see Section 6.1). The time taken for such a 
survey is just the time taken to image down to the median 
redshift for a given telescope, multiplied by the number of 
clusters. Note that the cumulative total number of haloes, 
A/"(> M,z m ), depends on the median redshift of the sur- 
vey, z m ; the upper scale on Figure 9 assumes z m = 0.7 and 
Zmax = 1.5. However, comparing with Figure 5 we see that 
varying z m has only a small effect in the number of haloes 
above M = 10 13 Mq, but does change the total number for 
masses below this. 

Figure 9 implies that by imaging only 60 of the most 
massive clusters (lO lj Af0) in a hemisphere to z m = 0.9 
(R = 24.7) in five bands and combining with the 4-year 
WMAP, one could reach an accuracy of Awo = 0.50, af- 
ter marginalizing over all other parameters, including w a - 
This seems a viable strategy, a factor of 2 improvement 
on 4- year WMAP given a marginalization over w a - If w a 
is fixed at w a = then the marginal error on wo reduces to 
Awo ~ 0.25, a factor of 2.5 improvement on 4-year WMAP, 
marginalizing over other parameters. To rapidly image each 
halo in five bands with an 8-metre class telescope with a 
0.025 square degree field-of-view, such as with SuprimeCam 
on the Subaru telescope (see Broadhurst et al., 2005, for 
the use of Subaru in a lensing analysis) would take 10 to 20 
nights. 

Beyond this accuracy, there are diminishing returns for 
a pointed survey from the geometric test. To reach an ac- 
curacy of Awo ~ 0.05, one would have to image around 10 7 
haloes, with the number of galaxies scaling roughly as 

N « io a35/A ™° . (84) 

For a targeted survey, this seems an unfeasible task. The 
weakness in this relation is due to the fact that we have 
ranked haloes by mass, and while the number of haloes is 
increasing the mass per cluster, and hence lensing signal, is 
falling. By the time we are imaging the < 1O 13 M0 haloes, 
the shear signal is so weak as to no longer contribute to a 
significant measurement of dark energy. 



© 2005 RAS, MNRAS 000, 1-30 



12 A.N. Taylor, T.D. Kitching, D.J. Bacon, A.F. Heavens 



Total Number (z m =0.7) 
60 10 5 5xl0 6 7xl0 7 



in 


V 


i 1 1 


i 

■ 


ci 




v 

A< \ 


z m =0.5 

7. =0.7 

m 


d 

T 1 






7 . =n.g 

m 


d 








in 
o 








d 









10 15 10 14 10 13 10 12 



Mass / M Q 

Figure 9. Variation of marginal error on wq with the mass of 
lensing cluster for a pointed survey with 20,000 square degrees 
accessible, and for z m = 0.5, 0.7 and 0.9. The dashed line has 
no S/N threshold, the solid line has a threshold condition set 
by equation (48). We assume a WMAP 4-year prior. Note that 
the cumulative total number of galaxies depends on the median 
redshift, z m , (see Figure 5). Here it is calculated for z m = 0.7, 
but has little effect. 

These curves scale with the survey median redshift 
roughly as 

Awo(>M,z m )&Aw (>M,z m = 0.7)(jjj\ , (85) 

where the increase in accuracy arises due to the increase in 
number of background galaxies reducing the shot-noise, and 
the increase in available clusters reducing clustering vari- 
ance. This approximation fails for the most massive clusters, 
where imaging deeper does not help as we are clustering- 
limited. 



5.2 Time-Limited Survey Mode 

In contrast to a targeted observation mode, one could also 
use a large survey telescope with a wide field-of-view to con- 
struct a general wide-field survey, and extract haloes from 
this for the shear ratio analysis. In this case it makes sense to 
restrict the amount of telescope time one can allocate to such 
a survey. In the next Section we discuss the optimization of 
such a survey. Here we shall assume the optimum survey pa- 
rameters and investigate how the signal is distributed across 
the mass spectrum of haloes. 

Figure 10 shows the cumulative gain in accuracy on wo 
as we add haloes of decreasing mass. We have marginal- 
ized over the remaining parameters (Q m ,Q v ,w a ), and cal- 
culated the Fisher matrix using the analysis of Section 3.3. 
We have assumed a fixed-time survey with median redshift 
of z m — 0.5, z m = 0.7 and z m = 0.9 (limiting magnitudes 
of r = 23, r = 23.8 and r = 24.7, respectively), combined 
with a 4- year WMAP prior (see Section 6.1). The lines for 
z m = 0.5 and z m = 0.7 cross at approximately 7 x 1O 14 M0 



this is interpreted as for a fixed time survey the optimal me- 
dian redshift varies slightly with the mass range of clusters 
used. As clusters of lower mass are included the optimal me- 
dian redshift behaviour converges so that z m = 0.7 yields the 
lowest error, note Figure 10 includes a 4-year WMAP prior. 
The area of each survey is of 38,400 square degrees, 10,000 
square degrees and 3,660 square degrees, respectively, ap- 
propriate for a survey with one, or more, 4-metre telescopes 
with a 2 square degree field-of-view (more than one would 
be needed for a z m = 0.5, 38,400 square degree survey). 
Note again that the upper scale (Total Number) for number 
of haloes depends on median redshift which is here assumed 
to be Zm = 0.7. The cumulative number of haloes is half 
of that for a given mass than for Figure 9 as the total area 
probed is half. 

Again we find that the largest haloes provide the largest 
contribution to the measurement of wo, with an error of 
Awo = 0.6 from the largest 30 haloes. The error has flat- 
tened off from 60 to 30 haloes. As with the targeted survey 
mode, the increase in accuracy for including smaller haloes 
has diminishing returns. However, given these haloes will 
already be in the survey, the limitation here is processing 
time, rather than telescope time. For a 10,000 square degree 
survey to z m = 0.7 (r = 23.8) we can reach an accuracy of 
Awo = 0.08 from the analysis of — 3x 10 6 haloes, down to 
haloes with M > 1O 13 M0. The majority of the signal (the 
steepest gradient in Figure 10) comes from the relatively 
numerous intermediate mass haloes with M ~ 5 x 10 14 Mq. 
Beyond this the signal-to-noise per cluster is too small to 
contribute to a measurement of too. 

For a time-limited survey, it is useful to parameterize 
how the uncertainty on wo scales with different telescopes 
and surveys by scaling the error with the fractional survey 
sky coverage, 

/sky = 40,000sq. deg.' (86) 
where A is the survey area, so that 

Aw = AwoC/sky = 0.25) ^ , (87) 

where 

Hence one can trade off telescope size and field-of-view (fov) 
with the survey time-limit, T, and the median depth, z m - 

To summarise Sections 5.1 and 5.2, while a reasonable 
sized pointed survey of around 60 of the largest clusters in a 
hemisphere combined with the 4-year WMAP results could 
rapidly measure two to around Awo = 0.50 in a short space 
of time, to improve the accuracy to a few percent would 
require an unfeasible amount of telescope time. However, a 
time-limited wide-field lensing and 5-band photometric red- 
shift survey could push the accuracy down to a few percent 
accuracy, for example Awo = 0.08 for a 10,000 square de- 
gree survey to z m — 0.7, with the analysis of the millions of 
medium sized clusters and groups (M > 1O 13 M0). 

Time-limited survey designs and their optimization in 
measuring w(z) are considered in further detail in Section 
6. 
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Figure 10. Variation of marginal error on wq with the mass of 
lensing cluster for a 10,000 square degree survey to z m = 0.7, and 
scaled to z m = 0.5 and z m = 0.9 with area of 10, 000(0. 7/z m ) 4 . 
The dashed line has no S/N threshold, the solid line has a thresh- 
old condition set by equation (48). We assume a WMAP 4-year 
prior. 



5.3 Area-limited Survey Mode 

A further distinct class of experiments, such as the LSST 
(see Tyson et al., 2002) and Pan-STARRS (PS4; Kaiser, 
2005), will repeatedly image an entire hemisphere (20,000 
square degrees) to a given median redshift; this is proposed 
to be done by stacking multiple images. In these cases the 
limiting factor is the amount of sky available to a given 
telescope, and time allowing for a given median redshift to 
be reached. Figure 9 shows that the marginal error on wo 
will vary as the median redshift of the survey as 



Aw = 0.07 



07 



/sky 

0.5 



-1/2 



(89) 



so that a z m = 0.9, A — 20, 000 square degree survey could 
image approximately 7 x 10 7 clusters between 10 12 and 10 15 
Mq, and achieve a marginal error of Awo = 0.05. A survey 
of this type is a viable alternative to the time-limited wide- 
field survey. 



6.1 Combining lensing with other dark energy 
experiments 

As well as gravitational lensing, there are other experiments 
which can probe dark energy, notably the CMB, Baryon 
Acoustic Oscillations (BAO) in the galaxy power spectrum, 
and the supernova Type la Hubble diagram. Individually 
each of these probes can probe dark energy, but suffer from 
degeneracies between wq and w a , and with other parame- 
ters. These degeneracies can be lifted by combining meth- 
ods. Since there are a number of different probes, these ex- 
periments can generate a number of combinations which can 
be compared for consistency and as a test for systematics. 
In addition, dark energy probes can be divided into meth- 
ods that probe just the geometric properties of the Uni- 
verse, and those that combine the evolution of mass cluster- 
ing and geometry. These may respond differently depend- 
ing on whether the apparent dark energy is vacuum energy, 
modelled as a fluid with negative equation of state, or a 
change in gravity on large scales. Again, with a combination 
of methods these possibilities can be explored. In this pa- 
per we shall only address the combination of methods under 
the assumption that the dark energy can be modelled by 
a negative-pressure equation of state. Finally, in this paper 
we do not consider the Integrated Sachs- Wolfe (ISW) effect 
directly, via cross-correlating galaxy surveys with the CMB, 
although this too can probe dark energy. 

The error analysis of a combination of independent ex- 
periments can simply be accounted for by summing over 
each Fisher matrix for a CMB, Type-la supernovae (SNIa) 
or a Baryon Acoustic Oscillation experiment (BAO) 



pTOT _ f gl tST-pP 



(90) 



where F p are the predicted Fisher matrices for each type of 
data. We examine three different dark energy probes, mo- 
tivated by experiments which will be contemporary with 
any experiment that could use the geometric test. The fidu- 
cial cosmological model used in the Fisher calculations for 
these CMB, BAO and SNIa experiments are: Q. m = 0.27, 
fi„ = 0.73, h = 0.71, cr 8 = 0.80, Q, b = 0.04, w = -1.0, 
W a = 0.0, the scalar spectral index n 3 = 1.0, optical depth 
to the surface of last scattering r = 0.09, the running of the 
spectral index, 



a,-, 



dn(k) 



(91) 



din A: 

with a n = 0.0, the tensor to scalar ratio r = T/S with 
r — 0.01 and the galaxy bias factor, b, which we set to 
b = 1.2. 



6 OPTIMIZATION FOR A WIDE-FIELD 
CLUSTER LENSING SURVEY 

Having investigated the source of the lensing signal which 
contributes to the measurement of Wo, and shown that a 
time- limited, wide-field survey can reach high-accuracy mea- 
surements of wo, we now proceed to optimize such a weak 
lensing and photometric redshift survey for a fixed time to 
measure the properties of dark energy from the geometric 
test. 



6.1.1 WMAP and Planck Surveyor CMB experiments 

Here we consider both a 4- year WMAP experiment and a 14- 
month Planck experiment, with predictions calculated using 
CMBfast (version 4.5.1, Seljak & Zaldarriaga, 1996). We 
have used a similar procedure to that outlined in Hu (2002) 
and Eisenstein et al. (1998). The Fisher matrix for a CMB 
experiment is: 

jpCMB V""* ST^ 9Cx l ,„ n-1 dCyi , nrA 

j =2^2^ -deT {CoYl)xY ^e- (92) 

l mi „ X,Y 
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darkCAM 



Area/sq degrees 


-'I'll 


^max 


^Bands 


10,000 


0.70 


1.5 


5 




Planck 


Band/GHz 


^bcam 


a T /10- a 


crp/10" 6 


44 


23' 


2.4 


3.4 


70 


14' 


3.6 


5.1 


143 


8.0' 


2.0 


3.7 


217 


5.5' 


1.3 


8.9 



WFMOS 



Area/sq deg re6S ^-bin fcmax 

/hMpc- 1 Bias 



2000 
300 



1.0 
1.0 



0.15 
0.15 



SNAP 



AW 



NsNIa. 



1.5 



17 



2000 



1.25 
1.25 



0.15 



Table 1. The main default values parameterising the Lensing, 
CMB, BAO and SNIa experiments considered in this paper. For 
further details of the surveys see Section 6.1 and Table 4. 



6.1.2 Combining with SNIa experiments 

We have calculated errors on parameters for SNIa exper- 
iments for the proposed SuperNova Acceleration Probe 
(SNAP; Aldering, 2005) supernovae experiment using a pre- 
scription similar to that outline in Ishak (2005) and Yeche 
et al. (2006). The Fisher matrix, defined by Tegmark et al. 
(1998) and Huterer & Turner (2001), is: 



F, 



SNIa \ ^ 



1 



dm(z) dm(z) 



^ [Am(z)] ! ddi 



09, 



(94) 



where 771(2:) is the apparent magnitude of a supernova at a 
given redshift and N z is the number of supernova bins in 
redshift. The apparent magnitude is related to the luminos- 
ity distance by 771(2) = M + 51ogioDz,(z) where Dl(z) = 
(Ho I c)(l+z)r(z) is the iJo-hidependent luminosity distance. 
The normalization parameter is Ml = M — 51ogio (Ho/c) + 
constant, where M is the absolute magnitude of a SNIa. 

The effective magnitude uncertainty in a given bin at 
a particular redshift, taking into account luminosity evolu- 
tion, gravitational lensing and dust and the effect of peculiar 
velocity uncertainty is given by (Kim et al., 2003) 



Am(z) = \ a?n + 



5<r„ 



cz In 10 



(95) 



where the scatter in peculiar velocities of cr„ = 500 kms 
is assumed, and the systematic limit 8 m — 0.02 (for a space 
based experiment). We use the standard set of 2000 sim- 
ulated SNAP supernova distributed in 16 redshift bins of 
width Az = 0.2 between redshifts 0.0 < z < 1.8 the number 
per bin taken to be the simulated sample from Yeche et al. 
(2006) and Virey et al. (2004). The full SNIa parameter set 
is (fi m , fi„, w , w a , h). 



where Cxi is the power for X = T,E, TE or B (Tempera- 
ture, E channel polarization, Temperature-E channel cross 
correlation and B channel polarization) in the I th multipole. 

The elements of the symmetric covariance matrix are 
given in Eisenstein et al. (1998). For example the TT element 
of the covariance matrix is given by: 



Cov ITT ~ 



where B\ 



(2£ + l)f aky 



(C T c + w T 1 B e 



(93) 



Gaussian beam window function B\ 



is a 

exp(-£(l + l)#b eam /81n2) and 6> beam is the full-width, half- 
maximum (FWHM) of the beam in radians. The inverse 
square of the detector noise level on a steradian patch for 
temperature and polarization is given by Wi = (#bcamO"0 _2 
where i = T,P. The sensitivity in piK per FWHM beam 
(AT /T or AP/T) is <Ti = 4 ix . 

For multiple channels the quantity vjB\ is replaced by 
the sum of this quantity for each channel. The values for 
#bcam and Oi for the various experiments were taken from Hu 
(2002) (Table I), the Planck parameters are shown in Table 
1. We have used a maximum £ m ax — 2000 and minimum 
imin = 10 in the summation over wavenumber. f 3 k y is set 
to 0.66 to simulate a typical galactic cut. 

The 11-parameter CMB cosmological parameter set is 
(fi m , Q v , h, <7 8 , ^b, i"o, w a , n s , r, a n , r = T/S). We do 
not include a marginalization over calibration of the CMB 
instrument. 



6.1.3 Combining with Baryon Acoustic Oscillations 
experiments 

We have modelled the errors on cosmological parameters 
for a BAO experiment, taking a WFMOS-type experiment, 
following Seo & Eisenstein (2003), Blake and Glazebrook 
(2003) and Wang (2006). The Fisher matrix for a BAO ex- 
periment can be approximated by 



pBAO = Y^[A\nP(k,z\ 



2 din P(k cB ,z) dlnP(k cS ,z) 
1 d6 t 06, 



(96) 



where P(k c g, z) is the linear matter power spectrum (see 
Eisenstein & Hu, 1998) at a redshift z including growth fac- 
tors for an arbitrary dark energy cosmology (see Linder, 
2003). The summation is over redshift bins, 2, and wavenum- 
ber k. k c g is an approximation to the observable wavenum- 
ber averaged over both radial and angular direction and is 
given by 



fcnff = k 



raa(z)H(2 



1/3 



(97) 



where the subscript fid refers to the comoving distance r(z) 
and Hubble parameter H(z) at the fiducial ACDM cosmol- 
ogy. The fractional uncertainty on the measurement of the 
power spectrum is given by 
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A\nP(k,z) = 2tt 



Vk' 2 Ak 



1 + 



nP(k,z) 



(98) 



where V is the volume of the survey. We assume nP = 1 for 
all surveys (see Seo & Eisenstein, 2003). 

The BAO survey assumed has two redshift slices centred 
on z = 1.0 (0.5 < z < 1.3) covering 2000 square degrees and 
z = 3.0 (2.3 < z < 3.5) covering 300 square degrees. The 
volume is calculated assuming the area and redshift ranges 
at the fiducial cosmology. 

We have also calculated the BAO prediction for a survey 
with an area of 10000 square degrees with a median redshift 
of z m = 0.7, using five redshifts bins with ranges centred 
upon z = 0.4 (0.3 < z < 0.5), z = 0.6 (0.5 < z < 0.7), 
z = 0.8 (0.7 < z < 0.9), z = 1.0 (0.9 < z < 1.1) and z = 1.2 
(1.1 < z < 1.3). To include the effect of photometric redshift 
uncertainty we add a radial damping term (see Zhan et al., 
2005) 



P(k eS ,z) -> P(fc cff , Z )e- c2CT ' (z)fc cV^ ld M 



(99) 



where cr~(z) is given by equation (58). 

Alternatively, in an effort to reduce the photometric 
redshift error, the matter distribution could be estimated 
by grouping galaxies into clusters each containing n perc i uster 
galaxies (Angulo et al., 2005). This would have the com- 
bined effects of decreasing the effective number density 
n — > n/n pcrc i ustor and decreasing the redshift error by aver- 
aging the error over the group a z (z) — > a z (z)/ ^n pclc i ustcl . 
We found for n pcrc i ust cr > 1 the marginal errors on wo and 
w a increase, since the effect of decreasing number density 
increases the fraction error on the power spectrum by more 
than the decrease in the photometric redshift error can com- 
pensate. Hence we find that using clusters for the BAO ex- 
periment here does not add to the results of the Planck CMB 
experiment. 

To ensure we are in the linear regime the maximum 
wavenumber used in all the surveys is k — 0.15 /iMpc -1 , and 
we use Ak = 5 x 10 -3 /iMpc -1 . The full parameter set used 
is (fl m , &v, h, bas, fib, wo, w a , n s , a„) where b is a bias factor 
parameterizing the mapping of the dark matter distribution 
to the galaxy distribution. An important assumption is that 
the bias is a constant on the scales probed. 



6.2 A Simplified Error Model 

Before considering the full problem of optimizing a weak 
lensing survey for the geometric dark energy test, it is use- 
ful to consider a simplified estimate of the parameter uncer- 
tainty, so that the more complex results can be understood 
in terms of simple relations between competing effects. The 
uncertainty on wq is roughly given by 



Aw 



wo jVNbN, 
where 

AT cl = A/fov 



_ / d\nR \ 
d V d In wq J 



N, 



1/2 



(100) 



(101) 



is the number of independent clusters or fields in the analy- 
sis, 



is the number of redshift bins behind the lens, where z m is 
here the median redshift of the survey and Az is the typical 
redshift error at that depth. The typical number of galaxies 
per bin is 



N, « fiN tot /N B N cU 



(103) 



where /; is the fraction of galaxies in the field behind the 
cluster, and N to t is the total number of galaxies in the sur- 
vey. The terms in this expression arise from two sources. The 
first, proportional to a e , is the intrinsic uncertainty per shear 
mode due to galaxy ellipticities, and can be beaten down by 
increasing the number of galaxies per redshift bin, or by av- 
eraging over more bins, or more clusters. The second term, 
proportional to C 77 is due to lensing by large-scale struc- 
ture in between the lens and the source bins, and can be 
reduced by increasing the number of redshift bins (with the 
approximation that each lensing bin is independent) and by 
averaging over independent clusters. The number of clusters 
in the sample scales with median survey redshift as 

..3.4 



N cl (M> 10 i4 M Q ) = 10^ 



(104) 



clusters per square degree, where we have cut the cluster 
sample off at 10 14 A4©, where we find the signal contributing 
to the measurement of wo vanishes (see Section 5). 

In general we will be interested in fixed-time surveys, 
where the survey time scales roughly as 



(105) 



where / s k y is the fraction of the sky covered by the sur- 
vey, Zm is the median redshift of the survey, and To is a 
time constant, the time to observe the whole sky to a me- 
dian redshift z m = 1 (i.e. to a limiting magnitude of 25 in 
the r-band; see equation 67), set by the telescope specifi- 
cations and number of observed bands. The time scales as 
the fourth power of the median redshift due to cosmological 
dimming effects and the need to detect the object against 
the sky background. As a concrete example we shall use the 
Canada- France-Hawaii Telescope (CFHT; Semboloni et al., 
2006; Tereno et al., 2004), which is a 3.6m telescope with a 
1 square degree field of view, integrating over 5 bands, for 
which To = 2 X 10 4 nights. We shall also assume a projected 
number density on the sky which scales with the median 
redshift of the sample as 



ri2(< z) — 30z m galaxies per square arcmin, 



(106) 



as measured from the COMBO-17 survey, an angle averaged 
shear-shear correlation function, 



C 77 = 10~ b z m b , 

and an intrinsic ellipticity dispersion 
a. = 0.3. 



(107) 



(108) 



With this simplified error model, we find the fractional error 
on wo scales as 

Aw 



= 0.062z~ 135 (l + 24.l4>Az 



,V2 



(109) 



N B w z m /Az 



(102) 



Wo 

The leading term here is due to shot-noise, while the second 
term in quadrature is due to large-scale sampling variance. 
Assuming we have ten redshift bins, so that Az = 0.07 is 
typical of the photometric redshift error, equation (109) min- 
imizes at Zm w 1.0. For a fixed-time survey we find that for 
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a shallow, low-z, wide area survey, the error on wo is domi- 
nated by shot-noise. Here the signal is not very large, and the 
number of background galaxies (and therefore combinations 
of background source planes) is too low. For a deep survey 
this becomes dominated by large-scale structure clustering. 
This occurs because we have to make the survey area smaller 
to compensate for the depth. Hence we have fewer clusters 
to average over and reduce the clustering noise. Both sources 
of noise increase with the size of the redshift error, Az. In 
the case of shot-noise this is again because we have fewer 
combinations of source planes to sum over. In the case of 
clustering noise-dominated there is a stronger effect because 
we have fewer source planes to average out the effects of 
clustering. 

6.3 Survey Optimization 

The optimizations discussed in the following Sections only 
include a CMB Planck experiment, the combination with 
further experiments is discussed in Section 7.2.3. For a weak 
lensing and photometric redshift survey on a given telescope 
for a set amount of observing time, the survey itself is char- 
acterized by the area, parameterized here by / s ky, the me- 
dian redshift, z m , of the survey in the band used for weak 
lensing (usually the r- or i-band) and the number of bands 
used for photometric redshift accuracy, -/Vbands- For a given 
number of bands we only have one free parameter, which we 
shall assume is the median redshift, z m . 

Our procedure is to vary z m , calculating the survey 
area by equation (78). With the galaxy number distribu- 
tion and number counts, we can calculate the Fisher matrix 
and hence the marginalized uncertainty on a measurement 
of u>o. Figure 11 shows the marginalized error on wo (as- 
suming a 14-month Planck experiment) for a D = 4m class 
telescope with a 2 square degree field of view for a variety of 
numbers of photometric bands. For example a 5-band sur- 
vey would be the case for, e.g., the Dark Energy Survey on 
the CTIO Blanco telescope or the darkCAM survey. The re- 
sults reflect our analysis of the simple analytic model. For 
a shallow, wide survey the lensing signal is not strong, the 
number of background galaxies is low and so the error on 
too is shot-noise dominated. The error on wo is poor beyond 
z m ~ 0.7, indicating that clustering noise is a strong effect. 
The small variation with the number of optical bands is due 
to the effect that, despite the marginal error of the geometric 
test decreasing, the intersection with the Planck experiment 
does not substantially change. This is investigated further 
in Section 6.4. 

The optimal survey is a, 5-band, 18,500 square degree 
survey with median redshift z m = 0.6, combined with a 14- 
month Planck survey. However note that the dependence 
on median redshift is shallow about the minimum and that 
the optimal survey when considering a figure of merit (see 
Section 6.6.2) is a 5-band, 10,000 square degree survey with 
median redshift z m = 0.7, so that from hereon, and in Sec- 
tion 7, we will use this as our fiducial survey design. 

6.4 Optical and Infrared surveys 

In the last few years multi-band surveys have started to open 
up the high redshift Universe. Hence it is now possible to 



in 




median 

Figure 11. The uncertainty on wo, marginalized over all other 
parameters, as a function of median redshift, z m , for a time- 
limited survey, assuming a prior from a 14-month Planck experi- 
ment. The survey area is A = 10, 000(,2 m /0.7) -4 square degrees. 
We set a lower limit of z = 0.5, which would correspond to a hemi- 
sphere. The solid line is for 5-band photometric redshift survey, 
the dashed line for 9-band and the dot-dashed line for 17-band. 
Note that the time constraint is only on the 5-bands, assuming 
that the other bands will come from other surveys. Note, see Sec- 
tion 4.4.1, that we take an upper redshift limit of z max = 1-5. 

combine 5-band optical surveys with 4-band infrared surveys 
for 9-band photometric redshifts. We can study the effect of 
varying the number of assumed additional bands available 
on the measurement of dark energy parameters by varying 
the photometric redshift error. Figure 12 shows the varia- 
tion of the accuracy on wo, marginalized over all the other 
parameters with a 14-month Planck experiment, as a func- 
tion of varying the accuracy of the photometric redshifts. 
We parameterize this by defining 

a z (z)=a (l + z). (110) 

A value of <tq — 0.05 is approximately appropriate for a 
5-band photometric redshift survey, while ao = 0.01 cor- 
responds to a 9-band (4-band infrared and 5-band optical) 
photometric redshift survey. For a 5-band survey (ao = 0.05) 
we find Awo = 0.075, while for a 9-band (4-band infrared 
and 5-band optical) photometric redshift survey (ao = 0.01) 
we find Awo = 0.071. Note this is distinct from a 9 band 
optical survey considered up until this point. 

If the photometric redshifts are degraded, for instance 
if fewer than five bands are available, the accuracy of wo 
is also degraded. By the time ao = 0.1 (for, say, 3-bands), 
the error has increased to Awo = 0.094. Note we have not 
included the effect of outliers here (see Section 8.4), which 
will degrade the signal further. 

We have found that using BAO to measure dark en- 
ergy from a photometric redshift survey is difficult as the 
damping term due to the photometric redshifts, effectively 
constraining the range of Fourier modes available to analyze, 
quickly reduces the amount of cosmological information that 
can be extracted. Figure 13 shows the variation of the error 
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Figure 12. The uncertainty on too, marginalized over all other 
parameters with a 14-month Planck experiment, as a function of 
photometric redshift accuracy, parameterized by a(z) = <ro(l + z). 
The normalization, ao scales roughly as the number of photomet- 
ric bands as an <x NT 1 , , where we find an = 0.05 for a 5- band 

u bands ' u 

photometric redshift survey and an = 0.01 for a 9-band (4-band 
infrared and 5-band optical) infrared and optical photometric red- 
shift survey. 



achievable using BAO from a photometric redshift survey, 
the error is simply the CMB error until <jo ~ 0.02 where the 
BAO constraint begins to improve the a 14-month Planck 
CMB error. To constrain dark energy using a photometric 
redshift survey many bands (possibly infrared) would be vi- 
tal over the whole redshift range to decrease the photomet- 
ric redshift error. As the redshift error becomes <J z {z) — > 0, 
as would effectively be the case for a spectroscopic survey, 
the geometric test constraints and the BAO constraints are 
comparable. 



6.5 Scaling results to other surveys 

To scale these results to other weak lensing surveys, equation 
should be used with a time calibration i.e. 



T 

Tit 



ZmO 



\D 



fov 
fov 



(111) 



The subscript refers to parameters time, median redshift 
and area of a survey on a telescope with certain diameter 
and field of view. The scaling applies between surveys with 
equal number of bands; for 5 bands the Canada-France- 
Hawaii Telescope Legacy Survey (CFHTLS) can be used, 
while for 17 bands COMBO-17 can be used. Although it 
can be naively assumed that the time for a given survey 
scales proportionally with the number of bands so that 
T () — > TnNi,n/Ni, where Nt, is the number of bands in the 
survey. 

One of two questions may arise. What is the error on wo 
(or w a ) that can be achieved given T nights on a given tele- 
scope, and freedom to choose the survey design? Or, given 
a survey of area A and median redshift z m what is the con- 



Figure 13. The uncertainty on wo, marginalized over all other 
parameters with a 14-month Planck experiment, as a function 
of photometric redshift accuracy, parameterized by an. The solid 
line are the geometric test constraints, the dashed line are the 
constraints using BAO from a 10, 000 square degree survey with 

Z m = 0.7. 





CFHT 


COMBO-17 


D(m) 


3.6 


2.2 


fov (sq deg.) 


3 


1 


N (bands) 


5 


17 




1.17 


0.7 


Area (sq. deg.) 


170 


1 


T (nights) 


500 


6 



Table 2. Default survey parameters for the 5-band CFHT Legacy 
Survey and the 17-band COMBO-17 survey. 



straint on too (or w a ) that can be achieved? Both of these 
questions can be answered using the information given here. 

If the field of view of the telescope is small enough so 
that only approximately one cluster will be observable per 
pointing then a targeting strategy should be used. In this 
case Figure 6 should be used so that given P pointings on a 
given telescope the appropriate marginal error can be pre- 
dicted. For a targeting strategy the time trade-off is deter- 
mined not by the total area covered but by the number 
of pointings. The number of pointings achievable given T 
nights to a redshift z m can be expressed, as 



P = 



T_\ (z m0 y /D_\ 
T ) V Zm J \D ) 



foVQ 

A ' 



(112) 



The achievable marginal errors from a targeting strategy are 
however limited due to the large amount of clusters which 
need to be observed for a tight dark energy constraint. 

Given the freedom to choose any wide-field surveys me- 
dian redshift, the optimal median redshift of z m ~ 0.7 is 
insensitive to the number of bands, when combined with a 
Planck prior (see Figure 10). Equation (111) should then be 
used, with the appropriate calibration, to calculate the area 
achievable given T nights. If the number of bands is 5, 9 
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Figure 14. The uncertainty on w(z), the dark energy equation of 
state measured at different redshifts, marginalized over all other 
parameters. For gravitational lensing combined with 14-month 
Planck experiment. This shows that the highest accuracy con- 
straint on w(z) occurs at z = 0.27 with Aw(z = 0.27) = 0.0298. 




0.5 1 1.5 

Z median 

Figure 15. The uncertainty on w(z), the dark energy equa- 
tion of state measured at different redshifts, marginalized over 
all other parameters for gravitational lensing combined with 14- 
month Planck experiment, and its dependence on median red- 
shift. The contours are lines of equal marginalized w(z) error, the 
numbers on the lines being the marginal error on that line. 



or 17 the appropriate line in Figure 10 then scales propor- 
tionally up (and down) with decreased (or increased) arial 
coverage from 10, 000 square degrees, for a 5 band survey 
i.e. Awq(A) = (0.075)(yl/10, 000)~\ If the number of bands 
is not shown in Figure 10 then Figure 11 can be used to 
find the minimum of the appropriate Awo vs. z m line (at 
z m = 0.7). This can then be scaled for a differing arial cov- 
erage as before. 

Given a fixed survey of area A and median redshift z m 
Figures 10 and 11 can be used in a similar way. Given the 
error in Figure 10 for a given median redshift Awo(z m ) 
the achievable error can be calculated using Awo(A) = 
Awo(zm)(W, 000/^4). In scaling between bands a similar in- 
terpolation between Figure 10 and Figure 11 can be per- 
formed. 



6.6 Constraining w(z) at higher redshifts 

6.6.1 Pivot redshifts 

As well as constraining the marginalized dark energy equa- 
tion of state, w(z), at z = (wo), we can combine the mea- 
sured accuracy of wo and w a to estimate the measured ac- 
curacy of w(z) at higher redshift. Here we can gain some in- 
formation by using the degeneracy between wq and w a (see 
Section 7), to find a redshift where the anti-correlation com- 
bines to minimize the error. Figure 14 shows the expected 
accuracy of w(z) — wo + w a z/(l + z) as a function of red- 
shift for a 5-band, 10,000 square degree survey with median 
redshift z m = 0.7, combined with a 14-month Planck sur- 
vey. The highest accuracy measurement occurs at z — 0.27, 
where w(z = 0.27) = 0.0298. This low-redshift pivot redshift 
for the geometric test is due to its insensitivity to w a . 

Figure 15 shows how the error on w(z) varies with both 
redshift, z, and with median redshift of the survey, z m , for 
the same time-limited survey. It can be seen that the mini- 



mization in the error at z ta 0.7 in Figure 11 is reproduced 
at the z = line (along the x-axis) of the plot, and Figure 14 
is reproduced by considering the variation in the error along 
the Zm = 0.7 line. It is clear that if one is concerned with 
optimizing a survey design to constrain the error on w(z) 
at an optimal redshift then there is little sensitivity to the 
survey design. This is due to the effect of intersection, that 
is even though the lensing only error may be varying the in- 
tersection of the lensing ellipse with the Planck experiment 
ellipse does not vary considerably in width (characterized by 
the width of the inner contour) or orientation (characterized 
by the value of z at which the error on w(z) minimizes). 

6. 6. 2 Figure of Merit 

A useful 'figure of merit' (Linder, 2003; Linder, 2006; Dark 
Energy Task Force, DETF 2006) in dark energy predictions 
can be constrained by considering the smallest area of pa- 
rameter space constrained by a given experiment. The dark 
energy equation of state can be written as: 



w(a) = Wi + w a (a,i — a) 



(113) 



where Wi = w(cn) and we have expanded around scale factor 
a;. The error on w(a) is: 

Aw(a) 2 = Awl + ( a i ~~ a ) 2 Awa + 2(dj — a)Cov(wi, w a )(114) 

where Cav(wi, w a ) is the covariance between Wi an w a (equal 
to the corresponding inverse Fisher matrix element). By tak- 
ing the derivative of this quantity the scale factor at which 
the error minimizes can be found 

Cov(li>i, W a ) 



— CLi + 



Awl 



(115) 



In the standard expansion in equation (6) a\ = 1 and the 
above expression reduces to the equation for the pivot red- 
shift. In this formalism the pivot redshift occurs when the 
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median 

Figure 16. The figure of merit as a function of median redshift, 
z m , for a time-limited survey, assuming a 14- month Planck prior. 
The survey area is A = 10, 000(z m /0.7)~ 4 square degrees. The 
solid line is for 5-band photometric redshift survey, the dashed 
line for 9-band and the dot-dashed line for 17-band. 

covariance between the Wi and w a is zero. This is equiva- 
lent to the pivot redshift in the formalism of equation (6). 
The ellipse at the pivot redshift is then the smallest ellipse 
constrained by a given experiment. Since this ellipse is de- 
correlated its area can be simply approximated by 

Aui(a p ivot) * Aw a . (116) 

This is the figure of merit used to quantify the performance 
of any given experiment: the smaller the figure of merit the 
tighter the constraints on the equation of state of dark en- 
ergy will be over a larger redshift range. Broadly it can vi- 
sualized by comparing Figure 15 and Figure 16, the figure 
of merit is minimized where the lowest contour in Figure 15 
is widest, this can be seen in Figure 16. It can be seen that 
the optimal experiment when considering the figure of merit 
is at a median redshift of z m = 0.7 for 5 bands. The figure 
of merit is shown for all considered experiments in Table 4. 



7 PARAMETER FORECASTS 

Having found the optimal survey strategy to measure the 
dark energy equation of state for a given experiment, we 
can now investigate the constraints on the full parameter 
space. Throughout we shall assume a 10,000 square degree 
5-band photometric redshift weak lensing survey with a me- 
dian redshift of z m = 0.7 (r = 23.8). 

In this Section we shall discuss dark energy parameter 
constraints from geometric lensing alone (Section 7.1), com- 
bined with the WMAP 4-year and 14-month Planck exper- 
iments (Section 7.2), and combined with a WFMOS BAO 
experiment and SNAP SNIa experiment in Section 7.2 and 
7.3. A table of the different surveys we have considered, and 
the predicted marginal errors on the dark energy parame- 
ters, is presented in Table 4. 



Using the full Fisher matrix formalism for parameters 
in a consistent cosmological model we can estimate the ac- 
curacy on a set of cosmological parameters for a given ex- 
periment, taking into account marginalization over all other 
parameters. The details of the Fisher analysis are discussed 
in Section 6.1. The 11-parameter cosmological parameter 
set we shall use is (fi m , f2„> h, as, Qb, wo, w a , n 3 , r, a„, 
r = T/S), with default values (0.27, 0.73, 0.71, 0.8, 0.04, 
-1.0,0.0, 1.0, 0.09, 0.0, 0.01). We shall compare and combine 
analysis with the results from a weak shear spectral analysis 
(e.g. Heavens et al., 2006) elsewhere. 

7.1 Parameter forecasts for the geometric lensing 
test alone 

On its own, the geometric test constrains a sheet in the 
likelihood space of (wo, W a , fiv, fim)- Figure 17 shows this 
plane in the 3-space of (wo, Q v , S2 m ), having marginalized 
over w a (light grey plane). The surface here encloses the 3- 
parameter, 1-er likelihood surface. The equation of this plane 
in the full 4-parameter space is 

X = 0.64to o - 0.31™ a - 0.350,, - 0.67O m . (117) 

For model parameters of wo = —1, w a = 0, fi m = 0.27, and 
f2„ = 0.73 this can be evaluated to give 

X = —1.08, (118) 

which can be measured with an expected accuracy of 

AX = 0.031. (119) 

If we fix wo = —1 and w a = 0, we can see that the geometric 
test constrains the degenerate line f2„ + 1.91fi m = 1.26. This 
can be compared with the CMB constraint on the density 
parameter plane of Q v + Q m = 1. 

We can project this onto a 2-parameter space, marginal- 
izing over all other parameters. Figure 18 show the 2- 
parameter, 1-er (68.3% confidence) likelihood contours for 
the parameter space of Q v , Q m , wo and w a - The lightest 
grey solid block is the constraint on parameters from the 
lensing geometric method only. Here again we see the large 
degeneracies between the geometric parameters. In particu- 
lar it is again clear that the geometric test is very insensitive 
to w a (see Section 2.3). The 1-parameter, 1-a marginalized 
parameter uncertainties can be found by projecting these 
contours onto each axis and dividing by 2.3. These are pre- 
sented in Table 4. 

7.2 Comparing and combining the geometric 
lensing and the CMB 

To lift the degeneracies in the geometric test we can combine 
our predictions with results expected from the CMB. Here 
we consider combining with the expected results from the 
4- year WMAP experiment. Below we shall compare with 
the results expected from a 14-month experiment with the 
Planck Surveyor. 

7.2.1 Combining with WMAP 

The parameter forecasts for a 4-year WMAP survey are 
compared and combined with the geometric test, allowing 
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Figure 17. 3D parameter space for a 10,000 square degree lensing 
survey to a median redshift of z m = 0.7 with 14-month Planck ex- 
periment, with no dark energy evolution. The volumes bounded 
by light and dark grey which represent the 1-cr parameter es- 
timations for weak lensing and a 14-month Planck experiment 
respectively. 

for spatial curvature, in Figure 18. The lightest grey ellipses 
are the geometric test alone, the darkest ellipses are the 
marginalized parameter forecasts for WMAP, while the cen- 
tral white ellipses, show the combined likelihood contours 
for the combined CMB and geometric methods. We have 
suppressed the amplitude of density perturbations param- 
eterized by erg, the Hubble parameter h, the optical depth 
t, and the tensor-to-scalar ratio r, which are also estimated 
by the CMB. We shall consider these parameters in Section 
7.2.2. 

Figure 18 illustrates the poor sensitivity of the CMB to 
wo and w a , but constrains the curvature of the model by the 
combination S2 m + Sl v . The response of the CMB to dark en- 
ergy comes mainly from the Integrated Sachs- Wolf (ISW) ef- 
fect. Combining the geometric lensing test and the CMB, we 
find the orthogonality of the two methods reduces the error 
on the dark energy parameters from Awo(WMAP) = 1.268, 
A™ a (WMAP) = 2.225 to Am (WMAP + GL) = 0.089 and 
An> a (WMAP + GL) = 0.714. There is also marginal im- 
provement in AQ m and Afi„ . The main improvement to the 
lensing analysis is the WMAP constraint on the curvature of 
the Universe in the Sl m , Hv parameter plane. To get a clearer 
picture of the orthogonality of the CMB 4-year WMAP and 
lensing geometric test results, we plot a 3-D view of the one- 
parameter, 1-cr parameter surfaces in Figure 17. This shows 
the wo, Sim, Civ parameter surfaces, marginalized over all 
other parameters, including w a . 



7.2.2 Combining with Planck Surveyor 

We can compare the information in Figure 18 from a 4-year 
WMAP experiment with that of a 14-month Planck Sur- 
veyor experiment, shown in Figure 19. While the Planck 
error ellipses (darkest grey) are considerably smaller than 
those of the 4-year WMAP, the degeneracy between wo 
and w a remains. On its own Planck can measure t«o to 
an accuracy of Awq — 0.502 and on w a to an accuracy 




Figure 18. Two-parameter, 1-cr (68.3% confidence) likelihood 
contours for geometric parameters for a 10,000 square degree lens- 
ing survey geometric analysis to a median depth z m = 0.7, com- 
pared and combined with the expected 4- year WMAP results. 



of Aw a = 1.86, with the main source of information from 
the Integrated Sachs- Wolfe (ISW) effect. Again the curva- 
ture of the model is well constrained by the CMB. Com- 
bining Planck with the geometric lensing test reduced the 
dark energy parameter uncertainties to Awo = 0.075 and 
Aw a = 0.326, a factor of ~ 7 improvement in the measure- 
ment of wo over Planck alone. 

The effect of the geometric tests constraints within an 
11-dimensional parameter space can be seen in Figure 20. 
All other parameters are marginalized over. Even though the 
geometric test does not place any direct constraint on the 
non-geometric parameters, we note that there is improve- 
ment in the normalization of matter perturbations, erg. This 
arises because erg, measured from the CMB is dependent 
on the parameters. Hu & Jain (2004) show the dependence 
of ag on other cosmological parameters, and in particular a 
constant value of w. In calculating the value of erg using dark 
energy dependent growth factors they find that the value of 
erg depends on a combination of dark energy parameters, 
they find an analytic expression in the special case of a flat 
Universe with constant w. These general arguments can be 
generalized to wo and w a using the growth factors given in 
Linder (Linder, 2003). An alternative parameter would be 
to use the horizon-scale amplitude of matter perturbations, 
5^, which is an independent parameter. We have chosen to 
use crs to compare with other analysis. The improvement on 
CMB parameters are summarized in Table 3. 
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Figure 20. The two parameter 1-<t (68.3% confidence) geometric constraints for a 10,000 square degree lensing survey to a median depth 
2 m = 0.7, with a 14-month Planck experiment in the 11-dimensional parameter space (Q m , Hu, h, ag, Of,, wo, w a , n s , T, a n , r). Other 
parameters for the CMB calculation arc marginalized over. 



7.2.3 Comparing and combining lensing with CMB, BAO 
and SNIa experiments 

Figure 21 shows comparisons between the geometric lens- 
ing, CMB, SNIa and BAO experiments for the geometric 
parameter set (fi m , fi„, wo, w a ). The broad, second lightest 
grey ellipses are for a SNAP-like SNIa experiment, the closed 
darker ellipses are for a WFMOS-like BAO experiment, the 
lightest grey is for the lensing geometric test, while the dark- 
est ellipse is for a 14-month Planck CMB experiment. The 
small white ellipse at the centre is the combined uncertainty. 
We have scaled the axes so that the full parameter degenera- 
cies can be seen. It is clear that allowing for spatial curvature 
and evolution of the dark energy opens up large degenera- 
cies in many of the experiments. Because of the large-data 
set and sensitivity of the CMB to parameters, the CMB pro- 
vides the strongest constraints alone. In particular we can 
see very similar degeneracies between experiments in the 
(Q v ,wo) plane, while there is some orthogonality between 
experiments in the (wo,w a ) plane. Combining experiments 
improves the constraints on all of the parameters. In par- 
ticular, allowing for spatial curvature we find Aioq = 0.043, 



and Auia = 0.108. We shall study the combination of exper- 
iments in more detail in Section 7.3. 

In Figure 22 we show the same set of parameters, 
but this time assuming spatial flatness. Again many of the 
largest degeneracies in each of the experiments remain. We 
also see clearly the insensitivity of the geometric lensing to 
w a , rendering it nicely orthogonal to the other experiments. 
Also note that the two geometric methods considered, the 
geometric test and the supernovae test constrain similar re- 
gions in the (wo, w a ) plane. Comparing Figure 21 and Figure 
22 it is clear that the assumption of flatness improves the 
marginal errors of the lensing, BAO and SNIa significantly, 
however since the CMB experiment constrain flatness to a 
high degree the overall combined constraints are broadly the 
same. This highlights the danger of assuming flatness, given 
that the marginal errors without a CMB experiment are 
drastically altered by this assumption. Given that some dark 
energy models involve variations to the Friedmann equation 
in non-flat geometries it is prudent to marginalize over spa- 
tially curved models. 

To illustrate further the orthogonality of the con- 
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Lensing 
Q CMB 
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Figure 19. Two-parameter, 1-cr (68.3% confidence) likelihood 
contours for geometric parameters for a 10,000 square degree lens- 
ing survey to a median depth z m = 0.7, with a 14-month Planck 
experiment. Note the change in the scale of the axes from Figure 
18, from hereon the remaining Figures will uses the scale intro- 
duced in this Figure. 



Parameter 


Planck only 


Combined 




0.0058 


0.0042 




0.0024 


0.0020 


h 


0.0088 


0.0070 




0.1002 


0.0383 


a b 


0.0011 


0.0008 


w 


0.5015 


0.0751 


W a 


1.8618 


0.3256 


n 3 


0.0034 


0.0034 


Ctn 


0.0062 


0.0056 


T 


0.0208 


0.0204 


r 


0.0079 


0.0077 



Table 3. Improvements on CMB Planck one parameter 1-cr, con- 
straints by adding the geometric test from a 10,000 square degree 
lensing survey to a median depth of z m = 0.7. 
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Figure 21. Two-parameter, 1-cr (68.3% confidence) likelihood 
contours for geometric parameters for a 10,000 square degree 
lensing survey to a median depth of z m = 0.7, combined with 
a 14-month Planck experiment, a WFMOS BAO experiment and 
a SNAP SNIa experiment. 1-parameter marginalized results arc 
tabulated in Table 4. 



straints from lensing, the CMB, BAO and SNIa, Figure 23 
shows a 3-dimensional plot of the likelihood contours in the 
(Q v ,wq, w a ) parameter space, marginalizing over all other 
parameters. We have plotted the 1-parameter, 1-sigma con- 
tours for clarity. 

7.3 Synergy of dark energy experiments 

It is interesting to compare the results of each of the dark 
energy experiments under the same conditions. In Figure 
24 we show the dark energy equation of state parameters 
(wo,w a ), marginalized over all other parameters including 
spatial curvature for each experiment in combination. 

Of all of the individual experiments considered the 
Planck CMB experiment on its own provides the strongest 
constraint on the (wo,w a ) plane, with the majority of the 
signal coming from the low-redshift ISW effect. However 
the marginalized uncertainties are still A^o = 0.502 and 
Aw a = 1.86. A SNAP-like SNIa experiment on its own pro- 
vides poor constraints in the wo, w a plane, due to the large 
degeneracy in Q m and Q v in models allowing curvature. 
This can be seen by comparing with Figure 22, but note 
that of the experiments considered, the supernova estimates 
are the only ones which include terms for extra systematic 
effects. Removing the extra systematic terms from the su- 
pernova estimates improves the constraints, when combined 
with a Planck CMB prior, by a factor of approximately 1.5 
to Aw = 0.094 and Aid„ = 0.318. A WFMOS-like BAO 
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Figure 22. Two-parameter, 1-a (68.3% confidence) likelihood 
contours for geometric parameters for a 10,000 square degree lens- 
ing survey to a median depth of z m = 0.7, with a 14-month Planck 
experiment, a WFMOS BAO experiment and a SNAP SNIa ex- 
periment, assuming spatial flatness with Q m + £2i> = 1. 




Figure 23. Likelihood contours in the 3-dimcnsional f2„, wo, Wa 
parameter space for geometric parameters for a 10,000 square 
degree lensing survey to a median depth of z m = 0.7, with a 
14-month Planck experiment, a WFMOS BAO experiment and a 
SNAP SNIa experiment, assuming spatial flatness, Q m + Q V = 1. 
1-paramcter, 1-sigma contours are used for clarity. 

experiment provides a narrow, but highly degenerate ellipse 
in the (wo,w a ) plane. This is due to the BAO experiment 
mainly constraining w(z) at the redshift of the nearest red- 
shift bin (in this case z = 1.0). Interestingly the BAO de- 
generacy is in a similar direction to the CMB degeneracy, 
presumably because a similar geometric effect is being mea- 
sured. Finally, the geometric lensing again has a large de- 




-2 2 

Figure 24. The combined marginal wo, w a constraints for two 
pairs of experiments. The experiments are a darkCAM lensing 
experiment and a CMB 14-month Planck experiment, a BAO 
WFMOS experiment and a SNIa SNAP experiment. Note that 
only the SNIa analysis contains terms for systematic effects. See 
Section 4.9 for details. 



generacy in this plane, but one which is different from the 
other experiments. The combination of pairs of experiments 
is very interesting. The combination of the geometric lensing 
and CMB puts very strong constraints on the dark energy 
equation of state and its evolution, reducing the uncertainty 
to Awq — 0.075 and Aw a = 0.326. Geometric lensing and 
SNIa yields Awo = 0.104 and Aw a = 0.699 while geometric 
lensing and BAO yields Aw = 0.128 and Aw a = 0.538. 
This provides us with three cross-checks with similar accu- 
racy. Looking at the dependency of each method, we see that 
both the geometric lensing, BAO and SNIa are all depen- 
dent on the geometry of the Universe, and so should give 
the same result, assuming that the wo,w a parameterization 
is valid. The CMB combines geometry with evolution of the 
potential field, particularly in the ISW effect. 

Looking at the other possible combinations without 
lensing we see there is a similar sensitivity to dark energy. 
We have already pointed out the degeneracy between the 
BAO and CMB constraints, and so their combination only 
marginally improves on the CMB alone. A similar result is 
found for combining CMB and SNIa. Finally BAO and SNIa 
provides an uncertainly similar to CMB alone. 

From this study, we conclude that the best pair com- 
binations come from combining geometric lensing with any 
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Figure 25. The marginal too, w a constraints for a combination 
of any three of the dark energy experiments. A darkCAM lensing 
experiment, CMB 14-month Planck experiment, BAO WFMOS 
experiment and a SNIa SNAP experiment. 
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Figure 26. The marginal uiq, w a constraints for a darkCAM 
lensing experiment and a CMB 14-month Planck experiment, a 
BAO WFMOS experiment and a SNIa SNAP experiment. 



of CMB, BAO or SNIa experiments, with Aw ~ 0.10 and 
Aw a ~ 0.50, and also the BAO and CMB combination. Mul- 
tiple combinations will also allow a degree of cross-checking 
for consistency. Other combinations are a factor of up to 5 
times poorer due to similar degeneracies between wo and w a - 
Combining three experiments, in Figure 25 we again see that 
strongest measurement of (wo,w a ) comes from combining 
the geometric lensing analysis with the CMB and BAO ex- 
periments with the uncertainty on wo and w a pushed down 
to Ami) = 0.047 and Aw a = 0.112. Adding the SNIa re- 
sults to this makes no a small difference (see Figure 26, but 
again recall that the SNIa is the only estimate to contain 
systematic effects). Again the all the three experiment com- 
binations all provide complimentary constraints in the (wo, 
w a ) plane. 



7.4 Complementary Figures of Merit and Pivot 
Redshifts 

The figure of merit and pivot redshift information can be 
represented as in Figure 27 so that both values can be seen 
simultaneously. The Figure shows a number of broad char- 
acteristics. As more experiments are added in combination 
both the pivot redshift converges to one mean value and 
the figure of merit decreases. The geometric test constraint 
forces the pivot redshift to lower values due to its unique de- 



I 9 

a O - 




pivot 

Figure 27. The figure of merit and pivot redshift for vari- 
ous experimental combinations. The combinations are labeled as 
L=Lensing, B=BAO, S=SNIa, C=CMB. Combinations of letters 
represent combinations of experiments. 



generacy whilst the BAO constraint forces the pivot redshift 
to higher values. It is also evident that the CMB constraint 
is not necessary for a low figure of merit (for example LBS). 
This Figure also shows how different combinations of surveys 
can probe dark energy at significantly different redshifts. For 
example the BC and LSC combinations both have a similar 
figure of merit with the BC combination z p i vo t = 0.76 and 
the LB combination z p i V ot = 0.28. 



7.5 The effect of changing the fiducial dark 
energy model 

The assumed fiducial cosmology has so far been a ACDM 
cosmology in which any derivatives in the Fisher matrix 
calculations have been about wo = —1.0 and w a = 0.0 
for the equation of state parameters. The effect of alter- 
ing this assumption is investigated here. We consider two 
alternative extremes which are just allowable by present 
constraints, dark energy models: a SUGRA (Super Gravity) 
model proposed by Weller & Albrecht (2002) represented by 
too = —0.8 and w a — +0.3; and a phantom model proposed 
by Caldwell et al. (2003) with w = -1.2 and w a = -0.3. 
To test the effect of changing our default dark energy model 
we re-run our Fisher analysis. 

Aswell as changing the point in parameter space about 
which the signal ratio is expanded in the Fisher matrix cal- 
culations the assumed fiducial dark energy model also affects 
the SIS to NFW scaling as a function of redshift and mass, 
as shown in Figure 6. It also affects the number density dis- 
tribution of haloes as a function of redshift and mass given 
by equation (49), when extending to arbitrary dark energy 
models we exchange the growth factor in equation (55) to 
the one given in Linder & Jenkins (2003). 

The effects of changing the default dark energy model 
are shown in Figure 28, where we plot the (wo,w a ) plane, 
fully marginalizing over all other parameters. The marginal 
errors for each experiment, and its degeneracy in parameter 
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Survey 


Area sqdeg 


z mcdian 


-^Bands 


Au>o 


Aui a 


z pivot 


Au>(z pivot ) 


Aui(z p i vot )Aui a 


Lensing 


darkCAM 


10000 


0.7 


5 


5.546 


31.132 


0.21 


0.972 


30.2471 


1 1 /~1 All i T~* 1 1 

darkCAM + Planck 


10000 


0.7 


5 


0.075 


0.326 


0.27 


0.030 


0.0097 


aarkUAM + UAU darkOAM 


10000 


0.7 


5 


0.459 


1.668 


0.01 


0.419 


0.6984 


darkCAM, 9 bands + Planck 


10000 


0.7 


9 


0.071 


0.311 


0.26 


0.029 


0.0089 


SNAP Lensing + SNIa + Planck 


1000 


1.38 


9 


0.073 


0.293 


0.31 


0.024 


0.0071 


All-Sky Space + Planck 


40000 


1.00 


9 


0.023 


0.146 


0.16 


0.012 


0.0017 


darkCAM+Planck+BAO+SNIa 


10000 


0.7 


5 


0.043 


0.108 


0.58 


0.018 


0.0019 


VST-KIDS+WMAP4 


1400 


0.6 


5 


0.227 


0.888 


0.19 


0.176 


0.1562 


CFHTLS(Wide)+WMAP4 


170 


1.17 


5 


0.282 


1.014 


0.29 


0.1663 


0.1687 


CMB 


4-year WMAP 








2.060 


3.612 


1.18 


0.758 


2.7379 


14-Month Planck 








0.501 


1.873 


0.367 


0.035 


0.0655 


BAO 


BAO WFMOS+Planck 


2000 


1.0 




0.070 


0.154 


0.78 


0.019 


0.0029 


SNIa 


SNIa SNAP+Planck 








0.142 


0.513 


0.37 


0.028 


0.0144 



Table 4. The table gives experimental parameters and marginalized cosmological parameter error forecasts for various surveys. 



space does indeed depend on the dark energy model. But the 
combined marginal errors do not significantly change. The 
main difference is manifest in the error on w a which increases 
for all methods as its value becomes more negative. This is 
simply due to the fact that a negative w a represents a dark 
energy scenario in which the dark energy density was less in 
the past (increasing in the future); so that the effect of dark 
energy on the expansion rate on observed galaxies (in the 
past) is less in these scenarios (and similarly the opposite 
effect for a positive w a )- 



8 ADDITIONAL SYSTEMATIC EFFECTS FOR 
LENSING 

In this Section we consider come potential systematic effects 
for weak lensing. To realize a ~ 1% measurement of u>o and 
w a from shear ratios not only do we need a large enough 
survey, but we must also be able to control systematics in 
a weak lensing survey to a high level. This will require con- 
trolling the systematics in the measurement of lens shear to 
A 7 ~ 1(T 5 . 



8.1 Image Shear Analysis 

The current generation of lensing surveys, with telescopes 
not specifically designed for lensing, induce 10% distortions, 
which can be corrected down to a net systematic of 0.01% 



(A7 ~ 10~ 4 ; Heymans et al., 2005). There has been exten- 
sive work into methods that can both diagnose and remove 
systematic errors from both intrinsic galaxy alignments (Hi- 
rata et al., 2004) and shear calibration errors. Mandelbaum 
et al. (2005) use a geometric test to diagnose systematic 
errors in the Sloan Digital Sky Survey (SDSS). 



8.2 Strong Lensing Effects 

There is a systematic effect in the strong lensing regime, 
where the reduced shear a, defined as 



(120) 

is measured from galaxy ellipticities. The ratio Rij now be- 
comes 



_ a _ 7i(l - Kj) 



(121) 



which, for the mildly non-linear lensing regime can be ap- 
proximated as 



Rii ~ — ~ (1 "4" K>i K j ) 

7j 



(122) 



Furthermore for a SIS k « 7. This was numerically tested, 
using observable clusters, and the amplitude of the correc- 
tion was found to be max(fti — Kj) ~ 0.15 and mean(Ki — Kj) ~ 
1.5 x 10 -3 . This numerical analysis implicitly assumes a ra- 
dius of 1 arcmin from equation 39, which is relatively narrow: 
tangential shear can be measured out to radii of at least 200 
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Is I 




-2 2 



Phantom 

w =-1.2 

w_=~0.3 



Aw : 0.05 
Aw :0.185 

a 



-2 2 



ACDM 
w =-1.0 
w = 0.0 



Aw n :0.043 
Aw :0.108 

a 



-2 2 



Quintessence 
w o =-0.8 
w> + 0.3 



Aw n :0.043 
Aw :0.099 

a 



Figure 28. The dependence of the marginal error on wo the 
assumed dark energy model, for a 10,000 square degree survey 
to a median depth z = 0.7, with a 14-year Planck prior, a BAO 
WFMOS prior and a SNIa SNAP prior. The errors quoted are 
marginal over all parameters. 



arcseconds see Gray et al. 2004. The largest source of this 
systematic error will be from the largest clusters, those pro- 
ducing the largest convergence, and as shown in Section 5 
the majority of the wo signal comes from clusters of inter- 
mediate mass for which we would expect this systematic to 
be smaller. 

Alternatively, one can construct a statistic which elim- 
inates the mass-dependence of 7/(1 — ft), such as the three- 
point statistic suggested by Gautret, Fort & Mellier (2000). 
This could be applied in the strong-lensing regime, again 
independent of the lens strength, and combined with the 
two-point geometric ratio test in the weaker lensing regime. 
We shall investigate this elsewhere. 



8.3 Cluster Substructure 

A further expected systematic is that arising from cluster 
sub-structure, which we assume is averaged over. The effect 
of including sub-structure can only increase the signal; as 
long as the mass map of a cluster can be accurately mea- 
sured the expected tangential shear signal can be modelled. 
One promising avenue which may yield information on sub- 
structure is flexion (see Bacon et al., 2005). In the low signal- 
to-noise regime (low galaxy number counts) in which a mass 
model may have to be assumed for a cosmological signal to 
be extracted then this systematic source of error will become 
important and the mass model will need to be accurately re- 
produced. However, in the high signal-to-noise regime where 
the number of available galaxies is such that ratio of the 
shears from the data can simply be taken this systematic 
source of error will not affect the analysis. 



8.4 The effect of photometric redshift outliers 

In any weak lensing photometric redshift survey there will 
be a sample of imaged galaxies that will not have photo- 
metric redshifts assigned. There will be several classes of 
objects, for some of which determining a photometric red- 
shift will be difficult. We test the effect of such 'outliers' 
here by assuming a population within the survey p2 that 
have photometric redshifts a p2 (z) = 0.5, that is they have 
practically no redshift information. 

There are two ways in which such a population can 
be used, either they are included in the sample of galax- 
ies used somehow, or they are discarded. If they are used 
then the galaxies can be treated as a seperate population, 
with of{z) = 0.5 and n p2 = (1 - A p2 )n where A p2 is the 
fraction of outliers in the total population, analysed inde- 
pendently and the constraints from the outliers added to 
the constraints from galaxies with good redshifts. Or, the 
effective redshift error distribution at a particular redshift 
2 can be modelled by the sum of two Gaussian distribu- 
tions, with errors a pl (z) and a p2 (z), the relative amplitudes 
of the Gaussians constrained so that A pl + A p2 = 1. Such 
a sum of Gaussians can be accurately modelled as an ef- 
fective Gaussian, see Blake & Bridle (2005), with an ef- 
fective width a e ff = ^AP^afiz)] 2 + AP 2 [af(z)} 2 . We in- 
vestigated varying the relative amplitudes of two Gaussians 
with a p (z) = a z (z), the original photometric redshift error 
from equation (58) and a second with a p2 (z) — 0.5. 

Figure 29 shows the effect of varying the fraction of to- 
tal population in oultiers, combined with a 14-month Planck 
prior, the solid line shows the constraints from treating the 
outliers as a seperate population analysed seperately, the 
dot-dashed line includes the outliers into a degraded popula- 
tion using an effective Gaussian. The second possibility, that 
of discarding the outlying sample, is investigated by simply 
reducing the surface number density by no — > (1 — A p2 )no, 
this is also shown in Figure 29 as the dashed line. 

All the methods for dealing with the outliers result in 
an increase in the marginal error on wo- As expected us- 
ing all the galaxies and treating the outliers as a seperate 
population has the smallest effect on the marginal error, by 
using the outliers the marginal error is less than when they 
discarded. The effective Gaussian method has the effect of 
decreasing the number of redshift bins in the survey that can 
be used while retaining the surface number density, thus de- 
creasing the signal. By discarding the outliers the number 
of redshift bins is retained while the surface number density 
is uniformly degraded. This shows that the signal is more 
dependent on the number of redshift bins, than the number 
density and that the strategy for dealing with outliers will 
be an important issue in future surveys. 



8.5 CMB Lensing 

When combining the shear ratio analysis with CMB mea- 
surements we have assumed that the weak lensing of the 
CMB by large scale structure and galaxy clusters can ef- 
fectively be ignored. Since the shear ratios do not contain 
any information about structure, there can be no correla- 
tion due to this. There may, however, be some correlation 
between lensing of the CMB and the noise term in the shear 
ratio method. We shall explore this elsewhere. 
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Figure 29. The dependence of the marginal error on wo on the 
amplitude of outliers with a trf 2 (z) = 0.5 for a 10,000 square de- 
gree survey to a median depth z = 0.7, with a 14- year Planck 
prior. The solid line combines the outliers constraint with the 
galaxies with good photometric redshift errors by adding the con- 
straints. The dot-dashed line shows the effect of combining the 
outliers with the rest of the population using an effective Gaus- 
sian. The dashed line discards the outlying sample of galaxies. 



9 CONCLUSIONS 

In this paper we have set out a new method for the analysis 
of the geometric shear ratio test for measuring the dark en- 
ergy equation of state, based on the measurement of shear 
ratios around individual galaxy groups and clusters. The 
shear ratio test is insensitive to the growth of structure, 
but sensitive to the geometry of the Universe, via the mat- 
ter and dark energy density and the dark energy equation 
of state. This approach allows one to apply the method 
to individual objects, rather than requiring the measure- 
ment of some other statistic such as the galaxy-shear cross- 
correlation function which may be noisy for small data-sets. 
The down-side is that the method is now contaminated by 
structure along the line of sight, which can be overcome by 
using many independent lines of sight. 

Of the parameters which govern the geometry of the 
Universe, or more properly the photon distance-redshift re- 
lation, the shear ratio is most sensitive to a constant dark 
energy equation of state, wo, and very insensitive to evolu- 
tion, parameterized here by w a - This can be understood as 
due to the shear ratios being sensitive only to the change 
in shape of the shear signal as a function of redshift. As 
w a parameterizes the high-redshift effect of the dark energy 
equation of state, its effects are "renormalized" away. This 
behavior is very different to other probes of dark energy, 
and so helps to break parameter degeneracies when com- 
bined with other probes. 

It must be emphasised that the Fisher matrix frame- 
work used in this paper may result in overly optimistic con- 
straints. Since the errors are calculated by expanding about 
a fiducial point in parameter space any higher order effects 
that may change the shape of the likelihood surface cannot 
be taken into account. The effect of varying the fiducial dark 



energy model, in Section 7.5, demonstrates that the errors 
are sensitive to the choice of the fiducial model. A concrete 
example of higher order likelihood effects can be seen in a 
3D cosmic shear analysis by comparing Fisher matrix calcu- 
lations of the (erg, flm) plane (for example in Heavens et al., 
2006) (predicting an ellipse) with the measured constraints 
from data (for example Kitching et al., 2006) which measure 
an extended curved constraint. These effects can be investi- 
gated by large simulations or by a Monte-Carlo type explo- 
ration of the likelihood surface, we leave such investigations 
for future work. 

To account for many of the sources of uncertainty in 
the method, we have developed a halo decomposition anal- 
ysis of the lensing dark matter distribution to model the 
signal from dark matter haloes over a range of mass scales 
and redshifts. We have also included the effects of shot-noise 
due to the random intrinsic orientation of each galaxy, pho- 
tometric redshift errors and the contribution of large-scale 
structure lensing to the error budget. We have also investi- 
gated in detail a model for the photometric redshift error, 
based on studies of the COMBO- 17 data-set, as a function 
of redshift, number of imaging bands and limiting magni- 
tude. The effect of a bias in the calibration and distribution 
of photometric redshifts with spectroscopic redshifts is also 
studied, and we find that we require some 10 4 galaxies with 
spectroscopic redshifts to control calibration issues. The lim- 
itations of observing the shear signal from the ground and 
space are also discussed, and we argue that without adap- 
tive optics ground-based lensing studies are seeing limited, 
suggesting that it will be difficult to use galaxies beyond 
z = 1.5. 

The halo decomposition analysis of the dark matter 
lenses has allowed us to probe the origin of the shear signal 
in different types of survey, taking a 4-metre telescope with 
a 2 square degree field of view as our default survey. These 
results can be scaled to any other telescope parameters. 

For targeted observations, where the time-limitation 
translates into the number of clusters and groups one can ob- 
serve to a given depth, we have shown that we only require 
around 60 of the largest clusters in a celestial hemisphere 
to constrain wq to around Atoq ~ 0.50, marginalizing over 
all other parameters, including w a , a factor of 3 improve- 
ment on 4- year WMAP given a marginalization over w a . To 
achieve a higher accuracy requires the imaging of an un- 
feasible number of haloes, and instead one should turn to 
a wide-field imaging and photometric redshift survey. We 
find for a 4-meter class telescope with a 2 degree field of 
view that with a 10,000 square degree, 5-band photometric 
redshift survey with median redshift z m = 0.7 (r = 23.8), 
we can expect to reach an accuracy of Awo ~ 0.07, again 
marginalizing over all other parameters including w a . Our 
results can be easily rescaled to other telescope types, and 
survey strategies. 

The halo decomposition allows us to deduce where the 
main signal comes from in both the targeted and survey- 
ing modes. In both cases a significant fraction of the sig- 
nal comes from the largest hundred clusters in each survey, 
reaching a sensitivity of Awo ~ 0.5, however the majority of 
the signal comes from the numerous (~ 10 5-6 ) M > 10 14 M© 
haloes which can push the accuracy up to Awo ~ 0.07. 

Having determined where the majority of the dark en- 
ergy signal will come from in a geometric shear ratio test, 
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we then investigate the optimization of such a survey, when 
combined with the expected results from the Planck Sur- 
veyor experiment. We find that for our fiducial telescope for 
a fixed-time survey, going shallower (z m < 0.7) over a wider 
area decreases the accuracy due to the drop in the number of 
available background sources and corresponding increase in 
shot-noise. Going deeper (z m > 0.7) over a smaller area in- 
creases the clustering noise, since we now have fewer clusters 
to average over. 

We have also studied the effect of varying the number of 
imaging bands to increase or decrease the photometric accu- 
racy. We find that when combined with Planck an increase 
from 5, 9 or 17 optical bands makes little difference to the 
optimal survey. The reason for this insensitivity to higher 
accuracy photometric redshifts is due to the integral nature 
of the lensing effect, and the weak effect when combined 
with another data-set. However decreasing the number of 
bands is expected to have a strong effect on the accuracy of 
the lensing survey as redshift information is lost. We discuss 
how our results can be scaled to other telescope classes and 
survey parameters. 

The dark energy parameters wo and w a can be com- 
bined to give an uncertainty on w(z) = wo + w a z/(l + z), at 
some optimal redshift. This combination helps distinguish 
where the survey is most sensitive to the dark energy equa- 
tion of state. In the case of our optimal lensing survey this 
is at z = 0.27 with Aw(z = 0.27) = 0.0298. Again, the 
reason for the low-redshift sensitivity to w(z) is due to the 
insensitivity of the shear ratio test to w a . 

Having optimized the lensing survey for the geometric 
test in combination with the expected results from the CMB, 
we have investigated the effect on the full set of cosmologi- 
cal parameters for the CMB and lensing. The geometric test 
constrains a narrow sheet in the (Q m , wq, w a ) parameter- 
space, which is nicely orthogonal to the CMB parameter con- 
straints. Here we show that the CMB mainly constrains the 
curvature of the model, while the geometric test constrains 
wo, and the combination constrain w a . 

We have also compared and combined the geometric 
shear ratio test with the expected results from an Baryon 
Acoustic Oscillation (BAO) experiment, such as proposed 
for WFMOS, and a supernova Type la survey, such as that 
proposed for SNAP. Here we have put all of the surveys 
(lensing, CMB, BAO and SNIa) on an equal footing, using 
the same curved background cosmology and the same dark 
energy model parameterization. We find that the degenera- 
cies in the geometric test, in particular the insensitivity to 
w a , are nicely orthogonal to all these other probes. Combin- 
ing the geometric test with the CMB, BAO or SNIa will yield 
accuracies of a Aiuo « 0.10 and Aw a m 0.5, and can be com- 
pared for systematics. An optimal combination is a geomet- 
ric lensing test, with the Planck CMB and WFMOS BAO 
experiment, yielding an expected accuracy of Atuo = 0.047 
and Aw a = 0.11. 

Finally we discuss some of the potential systematic ef- 
fects which could affect the predicted accuracy of lensing. 

In summary, the prospects of accurately measuring the 
dark energy equation of state and its evolution to high accu- 
racy over the next decade are very good. The key to this is 
the gravitational lensing geometric shear ratio test, which, 
due to its orthogonal degeneracies, can be optimally com- 
bined with a large range of other dark energy probes, such as 



the CMB, BAO or SNIa. In addition, gravitational lensing 
can also be a probe via two-point analysis, either from corre- 
lation functions or power spectra in redshift-space (Heavens 
et al. 2006 and Castro et al. 2005). Just as with lensing of 
the CMB, since the shear ratio analysis does not contain any 
information on structure, we can expect there to be little cor- 
relation between the two methods, even for the same survey. 
However, the shear ratio covariance may be correlated with 
the shear power. We shall explore combining these methods 
elsewhere. 
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APPENDIX A: COVARIANCE OF TANGENTIAL COSMIC SHEAR 

The cosmic tangential shear covariance averaged over a circular aperture is given by 

c^(9) = -1 / d 2 e fd 2 e' < 7 t(0)7t(0')> 



j d 2 6 j d 2 9' ( (71 (0)7i(0')) cos 2( P cos V + <72 W72 (0')>sin2<? sin 2tp') 

-i d 2 9 d 2 9' / -t—tx C" K (cos 2 2(^£ cos 2<^> cos 2ip + sin 2 2<y3£ sin 2ip sin 2ip')e^ ) 
A J J J 



dfl 
x £d£ 



1 { Jlajt-Q . 



Cg cos 2<^>f \ — j d 9e cos2ip 



cr 



2[1 - J (lfl)] _ Ji(^) ] 2 
f 2 (9 2 W J 



(123) 



where cosipe = l.6\ and 0i is the unit vector along one axis. 



APPENDIX B: BIAS IN ASSUMED PARAMETERS 

In this Appendix we show that for a Gaussian distributed likelihood function, the linear bias in a parameter, which we shall 
call 88i, due to a bias in a fixed model parameter (i.e., one whose value we have assumed and is not being measured), which 
we shall call Sxpj, is given by (e.g., Knox, Scoccimarro and Dodelson, 1998; Kim et al., 2004) 



(124) 
(125) 

(126) 



where F ee is the parameter Fisher matrix defined as 

Ft! = \?Ac- 1 d e l cc- 1 d]c + 2d%vC- 1 d%v T ) 

and F 9 ^' is a column matrix (for one bias parameter) defined as 

f m = ^(cr^ccr^c + dffic^df^ + df[iC- 1 d°v T ) 

which we will refer to as a pseudo-Fisher matrix between measured and assumed parameters. 

We begin with a Likelihood function, In L(8\ip), which depends on a set of free parameters to be determined by the data, 
0, and a set of fixed parameters which we assume are known, ip. If the are at their maximum likelihood values, 0o, then 

(0ilnL(0o|^)) =0 (127) 

where the derivative is in parameter space, and we have ensemble averaged over all possible data. 

We now ask what is the effect of displacing the fixed parameters. Expanding both ip and to first-order we find 

lnL(6»|V) =lnL(6> |Vo) +^ftlnL(6» jVo) +^^,,lnL(6> () |i/'o), ( 12 8) 
were d^,i is a derivative in the ^-parameter space. This displaced likelihood now maximizes when 

<cUnL(0|V)> = (a 4 lnL(0 o |V'o))+^(^lnL(0o|Vo)) +*&<&a*,ilnL(0 o |V'o)} = °- ( 129 ) 

We know that the unperturbed likelihood peaks at the maximum likelihood values, and by inspection we can see that the 
averaged second derivatives of the likelihood are the Fisher matrices. Hence we see 



89 j {d z d j \nL(0 o \tf> o )} = -^(diC^lnL^ojVo)) 

which with the definition of the Fisher matrices yields equation (124). 



(130) 
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